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Abstract. We study the local-in-time hydrodynamic limit of the relativistic 
Boltzmann equation using a Hilbert expansion. More specifically, we prove 
the existence of local solutions to the relativistic Boltzmann equation that are 
nearby the local relativistic Maxwellian constructed from a class of solutions to 
the relativistic Euler equations that includes a large subclass of near-constant, 
non- vacuum fluid states. In particular, for small Knudsen number, these solu- 
tions to the relativistic Boltzmann equation have dynamics that are effectively 
captured by corresponding solutions to the relativistic Euler equations. 
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1. Introduction and main results 

The special relativistic Boltzmann (rB from now on) equation provides a statis- 
tical description of a gas of relativistic particles that are interacting through binary 
collisions in Minkowski space, which we denote by M. The dynamic variable is the 
one-particle empirical measure F e > 0, which represents the average number of par- 
ticles of four-momentum P at each space-time point x G M. The four-momentum 
P of a particle of rest mass m is future-directecQ and satisfies the normalization 
condition P K P K = -m§c 2 , where the constant c denotes the speed of light. Conse- 
quently, we may view fasa function of time t G R, space iel 3 , and 3-momentum 
P= {P\P 2 ,P 3 ) GM 3 , witl 



\ 



A more geometric point of view is offered in Section |1.2[ where it is explained how 
to view F as function on the mass shell 

JTJt d = {( Xj P)eMx T X M | P K .P K = -mlc 2 , P is future-directed}, 

which is a submanifold of TM, the tangent bundle of M , and which is diffeomorphic 
to M 4 x R 3 . 

The relativistic Boltzmann equation (rB from now on) in the unknown F £ is 
(1.1) P K d K F e = ^C{F E ,F E ), 



where C(-, •) is the collision operator (defined in (1.10)), and the dimensionless 
parameter e is the Knudsen number. It is the ratio of the particle mean free 
path to a characteristic (physical) length scale. Intuitively, when e is small, the 
continuum approximation of fluid mechanics is expected to be valid. In this setting, 
we anticipate that the system of particles can be faithfully modeled through the 
use of macroscopic quantities, such as pressure, proper energy density, etc., whose 
evolution is prescribed by the equations of relativistic fluid mechanics, that is, the 
relativistic Euler equations. 

As a first rigorous step in this direction, we show that any sufficiently regular 
solution (n(t, x), 0(t, x), u(t, xy\ of the relativistic Euler system (rE from now on) 



satisfying the technical conditions (1.39) can be used to construct a corresponding 
family of classical solutions F e of the rB equation. Here, 9 denotes the the fluid 
temperature, n denotes the fluid proper number density, and u denotes the fluid 
four-velocity. Roughly speaking, the technical conditions are the assumption that 
9(t, x) is uniformly positive with only mild fluctuations, that n(t, x) is uniformly 
bounded from above and below away from 0, and that the spatial components of 
the four-velocity, namely u 1 (t, x), u 2 (t, x), and u 3 (t,x), are uniformly small. In 
Lemma |l.l[ we give a simple proof that these conditions are always satisfied near 
the constant fluid states. Under these assumptions, our estimates show that that as 
e — > + (in the hydrodynamic limit), the rB solution F e converges to the relativistic 



In the inertial coordinate system we use throughout this article, future-directed vectors P satisfy 
P° > 0. 

2 The formula for P° holds only in a coordinate system in which the spacetime metric g has the 
components g M „ = diag(— 1, 1, 1, 1). 
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Maxwelliar^\ M. associated to the solution of the rE system; see Section L4 for the 



definition of a Maxwellian, which is a function in equilibrium with the collision 
process; i.e., C(A4,A4) = 0. Thus, our results show that for small e, there are 
near-local equilibrium solutions to the rB equation whose underlying dynamics are 
effectively captured by the rE system. For the Newtonian Boltzmann equation, 
these are called "normal solutions" in Grad [24] . 

The main strategy of our proof is to perform a Hilbert expansion (see Section 



3.7) for F E . We write 



(1.2) F e = F + Ye k F k +e 6 F A 



k=l 



Re- 



inserting this expansion into ( 1.1 ) and equating like powers of e results in a hierarchy 
of equations. It turns out that Fq,...,F6, which do not depend on e, can be 
solved for: F must be Maxwellian, while Fi,...,F 6 solve linear equations with 
^homogeneities. Thus all of the difficult (and e— dependent) analysis is contained 



in the analysis of the remainder term Fr. e , which is carried out in Section 3.8 

Our methods and results can be viewed as an extension of the program initiated 
by Caflisch |8 , who proved analogous results for the non-relativistic Boltzmann 
equation and Euler equations [8]. We will utilize strategies from Caflisch as well as 



Guo [30] and Guo- Jang- Jiang [31] to perform the Hilbert expansion. We will also 

in 



use relativistic Boltzmann estimates from the work of the second author 45 47 
order to control the expansion. Additionally, we develop several necessary tools to 
study this problem in the setting of special relativity. In particular, we develop a 
mathematical theory for the kinetic equation of state, which is described just below. 
A more detailed discussion of the existing literature related to our result is located 
in Section |L8l 

Before proving the aforementioned results, we will sketch a proof of local exis- 
tence (in Section [2]) for the relativistic Euler equations. This local existence result 
ensures that there are in fact solutions to the rE system that can be used in the 
aforementioned construction. However, we mention upfront that during the course 
of our investigation, we ran into several technical difficulties that, to our surprise, 
seem to be unresolved in the literature. The first concerns the fundamental ques- 
tion of which fluid variables can be used as state-space variables in the rE system. 
In addition to the four-velocity u, there are five other fluid variables that play a 
role in the ensuing discussion: the aforementioned variables n and 8, together with 
the entropy per particle n, the pressure p, and the proper energy density p. In order 
to close the Euler equations, one must assume relations between the fluid variables. 



In this article, we assume that the three relations (1.27a I - (1.27c) hold between 
the five non-negative variables n,6,r),p, and p. As will be discussed below, these 
choices were not made arbitrarily, but are in fact satisfied by the macroscopic quan- 
tities n[A4], 0[yM], 77LW],p[A4], and p[M] corresponding to a relativistic Maxwellian 
Ai; these quantities are defined in Section [3~2] We emphasize that the relations 



(1.27a) - ( 1.27c[ ) are required in order for the rE system to arise from the 



rB equation in the hydrodynamic limit. Now it is commonly assumed that as 
a consequence of the three relations, any two of n, 8,r],p, p uniquely determine the 
remaining three. In particular, during our construction of the fluid solutions, we 
need to be able to go back and forth between the variables (n, 8) and the variables 

^These are also known as Jiittner distributions. 
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(jj,p), i.e., we need to be able invert the smooth map (n,z) — > (Sj(n, z), *P(n, z)) , 

2 

where z = "° g , ks > denotes Boltzmann's constant, and ,fj and ^ are defined in 
(1.29a I - (1.29b) below. However, we were unable to find a fully rigorous proof of 

Consequently, in Lemma |3.5| below, 



the invertibility of this map in the literature, 
we use asymptotic expansions for Bessel functions to rigorously verify the local 
invertibility of this map outside of a compact set. In particular, we show that 
the map is locally invertible whenever 9 is sufficiently large, and whenever 9 is 
sufficiently small and positive. Additionally, the numerical plot in Figure [I] (see 
Section [3~4| , which covers the compact set in question, strongly suggests that the 
map (n, z) — > (Sj(n, z), z)) — {i],p), is an auto-diffeomorphism of the region 
(0, oo) x (0, oo). This would imply that we can always smoothly transform back and 
forth between (n,9) and (rj,p) in the region of physical relevance, i.e., the region in 
which all of the quantities are positive; see Conjecture [l] in Section [2] below. Sim- 
ilarly, in Lemma 3.6 we rigorously prove that outside of the same compact set of 
temperature values, there exists a kinetic equation of state p = fkineticij], p), which 
gives the fluid pressure p as a function of the entropy per particle rj and the proper 
energy density p. 

A related issue is the fact that in order for the rE system to be hyperbolic 
and causaQ under a general equation of state p = f(r),p), it is necessary and 



sufficient to prove that < 



ill 



< 1. We explain why the positivity of 



is 



i] 



needed for our proof of local existence in Remark |2 . 1 1 below, while the mathematical 
connection between the inequality i 

r v 

being less than the speed of light is explained in e.g. 
kinetic equation of state p = fkineudv, p), — 



< 1 and the speed of sound propagation 
Now in the case of the 



dp 



43 



can be written as a function 



v 



of 9 alone. It is possible to write down a closed form expression for this latter 



quantity (see equation ( 3.32 ) ) , but since the formula is a rather complicated one 
involving ratios of Bessel functions, we have only analytically verified the inequality 
< 1 (again using asymptotic expansions for Bessel functions) outside 



< 



dp 



of the same compact set discussed in the previous paragraph; see Lemma |3.6| 
Therefore, the fully rigorous version of our local existence result is currently limited 
to initial data whose temperature avoids the compact set in question. However, we 
have numerically observed that in fact, the stronger inequality < 



f>h- 



dp 



< 



should hold for all 9 > 0; see Conjecture [2] in Section |2j and Figure [2] in Section 
3.5| This stronger inequality would imply that the speed of sound under the kinetic 
equation of state is never larger than \/l/3 times the speed of light. 

In view of these complications, when stating the hypotheses for our local ex- 
istence theorem (Theorem [l] in Section 1.7 1, we make careful assumptions on the 
fluid initial data that are designed to ensure that they fall within the regime of 
hyperbolicity, and within a regime in which the aforementioned map (n, 9) — > (rj,p) 
is invertible with smooth inverse. However, if our two conjectures are in fact cor- 



rect, then the relations (1.27a) - (1.27c) imply that many of these assumptions are 
automatically verified whenever the fluid variables are positive. Aside from these 
complications, our local existence theorem is a standard result. However, there are 
several additional aspects of it that are worthy of mention. First, we avoid the use 



By causal, we mean that the speed of sound is less than the speed of light. 
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of symmetrizing variables in our proof. We instead use the framework of energy 



currents, which was first applied by Christodoulou to the rE system in 11 , and 
which was later expounded upon by the first author in [ 43] . We also remark that 
our local existence result only applies to initial data with proper energy density 
p that is uniformly positive. In particular, we avoid addressing the complicated 
issue of the free-boundary problem for the relativistic Euler equations. A related 
comment is that our local existence result produces a spacetime slab [0, T] x M 3 on 
which the uniform positivity property is preserved. We remark that in view of the 
assumptions on n, 9 mentioned near the beginning of the article, we will only study 
fluid solutions belonging to compact subsets of the regions of interest to us, that 
is, regions where the maps and *}3 are rigorously known to be invertible. Thus, 
on such compact subsets, the uniform positivity of p is an automatic consequence 
of the continuity of the map (n, z) p, which is implicitly defined by the relations 



(1.27a I - (1.27b) 



1.1. Notation and conventions. We now summarize some notation and conven- 
tions that are used throughout the article. M denotes Minkowski space, while 93? 
denotes the mass shell. In general, Latin (spatial) indices a, b, j, k, etc., take on the 
values 1, 2, 3, while Greek indices k, A, p, u, etc., take on the values 0, 1, 2, 3. Indices 
are raised and lowered with the Minkowski metric g^„ and its inverse (g~ 1 Y v ■ For 
most of the article, we work in a fixed inertial coordinate system on M, in which 
case 

(1.3) ^ = Or 1 )"" = diag(-l, 1,1,1), 

and P K Q K = —P°Q° + $Z 0=1 P a Q a - Here and throughout, we use Einstein's sum- 
mation convention that repeated indices, with one "up" and one "down" are summed 
over. 

When differentiating with respect to state-space variables, we use the notation 



d> 



u\v, 



to mean partial differentiation with respect to the quantity U while V is held 
constant. We define x = {x\x 2 ,x 3 ), P = (P\P 2 ,P 3 ), P° = (m 2 c 2 + \P\ 2 )^ 2 , 
and |P| 2 = f X^i^ ) 2 ■ Furthermore Q,P',Q' are treated similarly. We also 
define P = f (P°) _1 P, and similarly for Q. We use the symbol 



^=(^'c^'^) =(9l ' 52 ' 53) ' 

to denote the spatial coordinate gradient. The Sobolev norm || • \\ h n of a Lebesgue 
measurable function f{x) on M 3 is defined in the usual way: 



1/2 



.\a\<N 



where d& = d^d^d^ 3 , a = (7^1,^2,^3) is a spatial coordinate-derivative multi- 
index, and \a\ — n\ + ni + n^. Here and throughout, we use the abbreviation 
h n dcf ff JV( K |)_ We also define the L 2_ inner product of tw0 f unc tions F(P),G(P) 
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as follows: 



(1.5) (F,G)p= [ F(P)G(P)dP. 



For brevity, we sometimes write (•, •) = (■, -)p. The L (Mp) norm is denoted | • I2. 

-2 

J x\P 



Similarly, we define the £ 2 .p inner product of two functions F(x, P), G(x, P) as 



(F,G) £; p= / f F(x,P)G(x,P)dxdP. 

Jr% Jm.% 

We denote the corresponding norm by ||/|| 2 = H/HiTO^R^ = ||/||l2(r| x r|)- We 
furthermore define the norm 

Halloo = ess sup xeK | !peM | \h(x,P)\. 

For each £ > 0, we also define the weight function we as 

(1.6) ^ = ^(P) d = f (l + |F| 2 ) £/2 . 

We then define a corresponding weighted L°° norm by 

WHoo.i = ess sup xeK | peR3p \w e (P)h(x,P)\. 

We define the norm of a Lebesgue measurable function f{x) over a measurable 
subset E C K| by 

\f\\Hg{E) = ( E \\ ds f\\Ll(E , 
K \S\<N 

\\f\\L UE) =[j E \f\ 2 dxy 2 , 

and similarly for the other norms and inner products over a subset. If X is a 
normed function space, then we use the notation C J ([0, T],X) to denote the set 
of j-times continuously differentiable maps from (0,T) into X that, together with 
their derivatives up to order j, extend continuously to [0, T]. We sometimes use the 
notation A < B to mean that there exists an inessential uniform constant C such 
that A < CB. Generally C will denote an inessential uniform constant whose value 
may change from line to line. For essential constants, we always write down their 
dependence explicitly. 

1.2. Lorentzian geometry and the mass shell 9Jt. In this article, we primarily 
work in a fixed inertial coordinate system on M, which is a global rectangular 
coordinate system {x M } M= o, 1,2,3 in which the spacetime metric has the form 



where 



(1.3). Note the sign convention of (1.3). This is the most common sign convention 
found in the relativity literature, but it is opposite of the sign convention that is 
sometimes found in the relativistic Boltzmann literature. 

Our coordinate system {x^} M= o,i,2,3 represents a special choice of a "space-time" 
splitting. We identify x° with ci, where c is the speed of light and t is time, while 
we identify (x 1 , x 2 , a; 3 ) = x with a "spatial coordinate": 

X — ^ X . X j X j X ^ — ct/ . X ^ . 
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We often work with the coordinate t rather than x°. Note that do = \dt. 
The mass shell is a subset of 

TM = U xeM T x M, 

the tangent bundle of M. In the following, we use an inertial coordinate system 
{x^, F"}^„ = o,i,2,3 on TM, where { 2^1^=0, 1,2, 3 is the inertial coordinate system on 
M, and 

t x m = ^y k ^ | (y\y 2 ,y 3 ,y 4 ) e m 4 J . 

In the above expression, (Y 1 , Y 2 , Y 3 , Y A ) are the coordinates of vectors in T X M ~ 
R 4 relative to the basis { J^-n I , J^r I , I , J^s 1 i . The mass shell 9Jt is defined 

L ox u \x ' ox L \x 1 ox z \x ' ox^ \x J 

to be 

M = {(x, P) e TM I P K P K = -rale 2 and P° > 0}. 
Let us also define 



m x = 371 n T X M, 



and 



0(P) d ^ f P° = J m 2 c 2 + \P\ 2 . 



It follows that the map $ : R 3 — > VJl x defined by 

$(P\ P 2 , P 3 ) d = f {^{P\P 2 ,P 3 ),P\P 2 ,P 3 ), 

is a diffeomorphism, and we can use it to put coordinates on 9Jl x ; i.e., any element of 
m xi viewed as a submanifold of T X M, has components (P°, P 1 , P 2 , P 3 ) = (4>(P),P) 
relative to our rectangular coordinate system. It follows that 

{x^, P- , } M =o,i,2,3j=i,2,3 is a global coordinate system on DJl. If (x, P) is an element 
of SOT, then we often slightly abuse notation by identifying (x, P) with (t, x,P). We 
similarly identify F(x,P) with F(t,x,P). 
We recall that there is a canonical measure 

(1.7) dflg ^ V^dctQ7)|dP, 

associated to <?, the first fundamental form of DJl x . We remark that g is Riemannian 
since 1R X is a spacelike hypersurface in T X M. This measure will allow us to define 
(in a geometrically invariant manner) integration over the surface dJl x . Recall that 
since $Jl x is (relative to the coordinate system {Y v } v= q\ 2 3 011 T X M) the level set 

SWx = {(P ^ 1 ^ 2 ^ 3 ) € T^M I P° - ^(P 1 , P 2 , P 3 ) = 0}, 

it follows that <? = where <I>*<7 is the pullback of g by $. Simple calculations 
imply that in our inertial coordinate system, we have 

gpJgpk gKX = ~Jpoy 



( L8 ) gj*= = ap ,- fl P fc g^ = -7^2 +(5 Jfc' (3,* = 1,2,3). 



Using (1.7) and (1.8), we compute that the canonical measure associated to g can 
be expressed as follows relative to the coordinate system P on dJl x : 

1 

P° l 



(1.9) d» 9 = -^dP, 



where we have used the fact that |det(g)| = (P ) -2 . We remark that (1.9) is valid 
only in an inertial coordinate system, and that integrals relative to the measure 
djig will play a central role in the definitions and analysis of Section 3.2 
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1.3. Hypotheses on the collision kernel. In order to state our hypotheses on 
the collision kernel, we introduce the following expression for the Boltzmann colli- 



sion operator, which is local in (t,x); we will elaborate upon it in Section 3.1 
(1.10) C(F, G) = P° [ v a{g, tf)[F(P')G(Q') - F{P)G(Q)]dQdu J , 



where we have suppressed the dependence of F and G on (t, x). Note that the 
collision operator acts only on the P variables. In the above expression, v = 
v (P,Q), the M0ller velocity, is defined by 



(1.11) 



v = v (P,Q) 



def C 



p 


Q 


2 i 


P Q 

P° X Q° 




Q° 


c 2 



C Qyfs 

4P°Q ' 



where x denotes the cross product inR . In ( 1.10), a is the differential cross-section, 
or the collision kernel. The relative momentum g{P, Q) is defined by 

def 



(1.12) g = ^{P- - Q-)(P K - Q K ) = V /-2(P«Q K + ?7i 2 c 2 ) > 0, 
while the scattering angle "$(P, Q, P' ,Q') is defined by 

(1.13) co S #^(P*-Q")(P' K -Q' K )/ Q 2 . 



Here the variables P' and Q' are defined in terms of P, Q below in (3.5 1. Finally, 
s(P, Q), which is defined by c 2 s = the energy in a center-of-momentum fram^j can 
be expressed as 

(1.14) s = -(P K + Q K )(P K + Q K ) = 2 (~P K Q K + m 2 c 2 ) > 0. 

Notice that s = g 2 + 4c 2 . We warn the reader that this notation, which is used in 



12 , may differ from other authors notation by a constant factor. Furthermore, u) 
is an element of S 2 (viewed as a submanifold of M 3 ), which can (with the exception 
of the north pole) be parameterized by the angles (t?, (p) £ (0, tt] x (0, 2ir], where i? 
is from above, and (p is an azimuthal angle. Relative to these coordinates, we have 
that du = sm-ddiSdip. 

The function a depends on the chosen model of particle interaction. For the 
remainder of the article, we assume the following: 

Hypotheses on the collision kernel: We assume that there are constants C\, 
C<i, a, andb such that the differential cross-section a, which is listed above in (1.10) 
and below in (3.4), satisfies the inequalities 

<t(q,0) < (d g a + C 2 g- b ) 0-0(1?), 

where 

• C\ and C2 are non-negative 

• There exists a number 7 > —2 such that < 00(1?) < sin 7, # 

• < a < min(2, 2 + 7), and < b < min(4, 4 + 7). 

We also assume that there exist constants C > 1 and (3 £ (—4, 2) such that 

(1.15) U p y /2 <v(p)<c(p°f /2 . 

center-of-momentum frame is a Lorentz frame in which + = P f f* + Q 1 ^ = (y/s, 0, 0, 0) 
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Here v(P), the collision frequency, is defined by 



def 



dQ I duj v cr(e,i?) J(Q), 
Is 2 



and J{Q) is the global relativistic Maxwellian defined by 
(1.16) J(Q) ^e-^/fte^), 



and 9m > is a constant. The lower bound in (1.15) is satisfied if there is a suit 



able lower bound for a. See 16 18 45 48 , and we refer specifically to 45 for more 
details. 



The upper bound in (1.15) is sufficient to deduce some of the estimates that we 



use below, such as Lemma 3.12 Lemma |3.13[ Lemma 3.14 and Lemma 3.15 Our 



hypothesis originates from the general physical assumption introduced in 16 ; see 
also [17] for further discussions. Standard references in relativistic Kinetic theory 
include (9j[l2j[l8j|44|[52j . 

1.4. Maxwellians. We now introduce the Maxwellians, a special class of functions 
on 9Jt that play a fundamental role in connecting the relativistic Boltzmann equation 
to the relativistic Euler equations. Given any functions n = n(t,x), 9 = Q(t,x), 
u p. _ u ^{t^x) on M such that n > 0, 9 > 0, u° > 0,u K u K = — c 2 , we define the 



corresponding Maxwellian M = M.(n, 9, u; P) as follows: 



(1.17) 



M = M(n,9,u;P) 



AnniQC 3 K 2 (z] 
where the dimensionless variable z is defined by 



exp 



(1.18) 



dcf m c 



ks is again Boltzmann's constant, and K%(z) is the Bessel function denned in (3.43) 
It can be shown (see e.g. [l2j Chapter 2]) that 



(1.19) C(F,F) = 



F is a Maxwellian of the form J\4(n, 9, u; P). 



For this reason, a Maxwellian M is said to be in local equilibrium with the collision 
process. If n(t, x), 9(t, x), and u^it, x) are constant valued, then M is said to be in 
global equilibrium. Note also that the function J(P) appearing in (1.16) is a global 
Maxwellian. It will play a distinguished role in the analysis of Section |3.8| In 
fact, we will assume that the local Maxwellian corresponding to the fluid solution 



is uniformly comparable to powers of J(P); see (1.39) 



Remark 1.1. Although C(Ai,Ai) = 0, it is not true in general that M. is a solution 
to the rB equation ( |1.1[ ); i.e., in general, P K d K (M.{n, 9, u\ P)) ^ 0. Nevertheless, 
our main result (Theorem [2]) shows that under certain assumptions, including that 
(n, 9, u) are solutions to the rE system, there is a solution of ( 1.1 ) near M.. 



1.5. The energy-momentum tensor and the particle current for rB. We 

now define TboUz\F\, which is the energy-momentum-stress-density tensor (energy- 
momentum tensor for short) for the relativistic Boltzmann equation, and IboUz[F], 
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which is the particle current. Given any function F(P), these quantities are defined 
as follows: 



(1-20) 



f - dP 

TZu z \F]=c P^F(P) — , (0<^<3), 
Jml P 



''<•/ (0<M<3). 



Whenever there is no possibility of confusion, we abbreviate TboUz = Tb Uz[F], 
and similarly for the additional quantities depending on F that appear below. 



The conservation laws can be summarized as follows (see Lemma 3.2): whenever 
F(t, x, P) is a classical solution to the relativistic Boltzmann equation, the following 
conservation laws hold 

fl21 , d K (TZJ F })=^ («<M<3), 

[ ■ ' d K {I Boltz [F])=Q. 

It is explained in the next section that whenever F = A4(n,9,u; P), and (n,9,u) 
are a solution to the rE system, then the above conservation laws for F hold even 
though F need not be a solution to the rB equation. This fact will play an important 
role during our discussion of the relationship of the rE system to the rB equation. 

1.6. The relativistic Euler equations and their relationship to the rela- 
tivistic Boltzmann equation. In this section, we recall some basic facts about 
the rE equations in Minkowski space. This is intended to serve as background for 
Theorem [T] which is stated in Section [L7j and proved in Section [2j For a detailed 
discussion of the rE system, we refer the reader to Christodoulou's survey article 



10 ; here we only provide a brief introduction. Our other goal in this section is to 
illustrate some of the formal correspondences between the rE system and the rB 
equation. These provide a heuristic basis for the expectation that the rE system 
should emerge from the rB equation in the hydrodynamic limit. 

The rE system models the evolution of a perfect fluid evolving in a Lorentzian 
spacetime. In Minkowski space, the spacetime of special relativity, they are 

(1.22) 



9*^ = 0, (0< M <3), 



Here the energy-momentum-stress-density tensor (energy-momentum tensor for 
short) for a perfect fluid has components 

(!-23) T^ uid A = f c- 2 (p + p)u^ +p(g- 1 r, (p, v = 0,1, 2, 3), 

where p > is the proper energy density, p > is the pressure, and u is the four- 
velocity. The four- velocity is a future-directed (i.e., u° > in our inertial coordinate 
system) vectorfield that satisfies the normalization condition 

(1.24) u K u K = -c 2 . 

The vectorfield Pfi uid is the particle current, it is proportional to the four-velocity: 



^„ id = n< ( M = 0,l,2,3). 

The quantity n > is the proper number density. All of these quantities are 
functions of (t, x) e K x R 3 . 
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By projecting the first divergence in ( 1.22 1 in the direction parallel to u and onto 
the g— orthogonal complement of it, we can rewrite (1.221 (omitting some standard 
calculations) in the well-known form 



(1.25a) 
(1.25b) 
(1.25c) 



d K (nu K ) = 0, 
u K d K p + (p + p)d K u K = 0, 
(p + pKcW +U^d K p = 0, 



(^ = 0,1,2,3), 



where II, the two-tensor that projects onto the g— orthogonal complement of u (i.e. 
U tJ,K u K = 0) has the components 



(1.26) 



IP 



— c u h u 



+ GT 1 ) 



(M^ = 0,l,2,3). 



The above equations are redundant in the following sense: if (1.24) and (1.25b) 
hold, and if (1.25c) holds for p = 1,2,3, then it follows that (1.25c) also holds 
when p = 0. 

The equations ( 1.22 ) are not closed because there are more unknowns than equa- 
tions. In order to close the rE system in a manner compatible with the rB equation, 
we will make use of the additional fluid variable 77, a non-negative quantity known 
as the entropy per particle. We now make the following assumptions, which are of 
crucial importance: we assume that the fluid variables n, 9, 77, p, p are bound by 
the relations 

p = ksn8 = rrioc 2 — , 



(1.27a) 
(1.27b) 

(1.27c) 



2 Ki(*) 
K 2 (z) 



n = 47re 4 mQ c 3 /i 3 exp 



3p, 



k B 



cxp 



V K,(z)J 



K 2 {z), 

where kg > is Boltzmann's constant, h > is Plank's constant, and z is defined 
in (1.18). Furthermore, the Kj(-) are modified second order Bessel functions, which 
are defined in Lemma 3.7 The origin of these relations, which are fundamentally 
connected to the properties of Maxwellians, is explained later in this section. 

Using the above relations, we can deduce the local solvability of any one of the 
variables n, 8, 77, p, p in terms of any two of the others whenever we know that 
the necessary partial derivatives are non-zero (knowing this allows us to apply the 

implicit function theorem). In particular, whenever > 0, we can locally solve 
for p as a function j 'kinetic of 77 and p : 

(1-28) p = fkineticiViP)- 

We refer to f kinetic as the kinetic equation of state. Note that this equation of state 
is discussed in Synge [52] . For future use, we denote by Sj and *P respectively the 
smooth maps from (n,z) to 77, p induced by the above relations: 



(1.29a) 
(1.29b) 



j] = Sj(n, z) = ks In < (Aire tuqC h 



.- 3 „-itfa(*) 



cxp I z 



KM 
K 2 {z) 



p = y{n,z) 



moc 2 nz . 



In the next section, we sketch a standard proof, which is based on energy esti- 
mates, of Thcorcm[T] i.e., for local existence for the rE system under the relations 
(1.27a) - (1.27c). During the proof, we work with the unknowns (77, p, u 1 , u 2 , u 3 ), 
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the reason being that a framework for deriving energy estimates in these variables 
via the method of energy currents has been developed; this is explained in detail 
during the proof of the theorem. To derive the energy estimates, we of course 
need an equivalent (for C 1 solutions) formulation the rE system (1.22) in terms of 
(tj,p, u , u 2 , u 3 ). To prove the equivalence of the systems, one needs the following 



identity, which is shown to be a consequence of the relations (1.27a) 
Proposition |3.4| below: 



(1.27c) in 



(1.30a) 



P + P 



dp 

dn 



Using ( 1.25a), ( 1.27a) - ( 1.27c I, and several applications of the chain rule, we deduce 
the following well-known version of the rE system: 



(1.31a) 
(1.31b) 
(1.31c) 

(1.31d) 
(1.31e) 



u k 8 kV = 0, 
u K d K p + qd K u R = 0, 
(p + p)u K d K u» + W* K d K p = 0, 



IP" = <rVu" + (g' 1 )^ , 



( M = 0,l,2,3), 
{H,v = 0,1, 2, 3), 



q = c 2 (p + p) 



dp 
dp 



Although our local existence theorem is proved using the variables (?/,p, u , u 2 , u 3 ), 
in order to construct the relativistic Maxwellian (1.17), which plays a fundamental 
role in our analysis of the rB equation, we require the availability of the variables 
(n, 9, u , u 2 , u 3 ). Thus, we need to be able to solve for (n, z) as a smooth function 
of (r],p). Remarkably, we could not find such a result in the literature. Thus, in 
Lemma 3.4 we rigorously show that we can locally solve for (n, z) in terms of (r),p) 
if < z < 1/10 or z > 70. Furthermore, based on numerical observations in the 
region 1/10 < z < 70, we make the following conjecture: 

Conjecture 1. The map (n,z) ($j(n, z), ?P(n, z)) is auto-diffeomorphism of the 



region (0, oo) x (0,oo), where the maps S) and are defined in (1.29a) - (1.29b) 



This conjecture is based on a numerical plot (see Figure 



0' 



in which z 5 

az 



appears to be negative for all z > 0. Examining the proof of Lemma pT4| it is clear 
that the negativity of would imply that the conjecture is true. We remark 



that by (3.24), analytically verifying the negativity of this quantity is equivalent to 
demonstrating the following inequality for all z > : 



(1.32) 



K 2 (z) \K 2 {z) 



\K 9 (z)J 



4 <o. 

z 



Another fundamental quantity in the system (1.31a) - (1.31e) is the speed of 



sound, the square of which is defined to be c 2 



It is a fundamental thermo- 



dynamic assumption that < 



< 1 for physically relevant equations of state. 



is required in order for 



As we alluded to in the introduction, the positivity of -rf- 

the rE system to be hyperbolic, while the upper bound of 1 implies that the speed 
of sound propagation is less than the speed of light. In Remark |2.1| we explain 
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exactly how we use the positivity in our proof of local existence, while the analytic 
and geometric connections between the upper bound of 1 and the speed of sound 
propagation being less than the speed of light is explained in e.g. [43]. Remarkably, 
there seems to be no rigorous proof in the literature that these inequalities hold 



for the kinetic equation of state (1.28) in every regime. Consequently, in Lemma 



3.6 we analytically verify that the equation of state ( 1.28 ) exists and satisfies these 



inequalities in the same regime discussed above, namely for < z < 1/10 and 
z > 70. Moreover, we conjecture that the following stronger statement is true. 



Conjecture 2. Under the relations (1.27a) - ( 1.27c[ ), p can be written as a smooth, 
positive function ofr],p on the domain (0,oo) x (0,oo), i.e., the equation of state 



(1.28) is well-defined for all (r),p) € (0, oo) x (0, oo). Furthermore, on (0, oo) x 

def dfkinetic(V: P) 



(0, oo), we have that 
(1.33) 



< 



dp 



dp 



< 



l 



Our conjecture is based upon a plot of {d z \ r) p)/{d z \ v p) (which can be expressed 

alone) that covers the region in question, i.e., the 



in terms of a function of z 



region 1/10 < z < 70. Our plot is labeled as Figure [2] of Section 



we note that according to equation (3.32 ), proving the inequality ( 1.33 ) is equivalent 



3.5 Furthermore, 



to proving that the following inequality holds for z > : 



(1.34) 



3 < 3 + z 



Ki(z) 
K 2 {z) 



'K 2 (2 



+ z 



\K2(Z)) 



o Ki(z) , „ 
6 K 2 (z) + z 



\K 2 (z)) 



< OO. 



is the middle term above ex- 



In this inequality (d z \ v p) / (d z \ n p) = ( g£y<e j 

pressed as a function of z. 

Now that we have provided a thorough discussion of the subtleties, we state our: 
Running Assumption: For the remainder of the article, we re- 
strict our attention to regimes in which any three of the five macro- 
scopic variables n, 9, 77, p, p can be written as smooth functions of 

the remaining two, and in which < I 2 (Vip)- To avoid bur- 

p v 

dening the paper, we do not point out this assumption 
every time it is made. Thus, under these running assumptions, 
any two of the five variables, together with the spatial components 
u , u , u of the four- velocity, may be considered to be state-space 
variables (i.e., they completely determine state of the system), and 
the phrase "a solution to the rE system" is a well-defined con- 
cept, independent of the choice of state-space variables. The role 



of Lemmas 3.5 and 3.6 then, is to guarantee that there is a large 
regime in which these conditions hold. Furthermore if Conjectures 
[l] and [2] are correct, then due to the relations (1.27a) - (1.27c), 
these conditions arc automatically verified. 
We now explain the connection between solutions to the rE system and solutions 
to the rB equation. First, given any relativistic Maxwellian A4(n, 9, u; P) as defined 



in (1.17), we can define an associated proper number density n[.M], temperature 
0[A4], entropy per particle rj[M], pressure p[.M], and proper energy density p[.M]. 
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The first two quantities are denned to be the first two arguments in M(n, 9, u; P), 
while precise definitions of the remaining three (as well as equivalent definitions of 



the first two) are given in Section 3.2 Furthermore, it can be shown that 



(1.35) 



1 fluid' 



fluid' 



Equation (1.35) should be interpreted as follows: the tensor obtained by substitut- 
ing p = p[M], p = p[M], and = I^[M]/n[M] into the definition (1.23) agrees 
with the tensor calculated using ( 1.20 ). The key point is the following: as shown in 



Proposition 3.3 the macroscopic quantities n[.A/J], 0[yM], ^[M], p[A4], p[M) satisfy 
the relations (1.27a) - (1.27c). Therefore we must use the same relations to close 



the rE equations 
Alternatively 
verify the rE equations (|1.31a) 



Alternatively if we are given smooth functions (77, p, it 1 , it 2 



3 ) of (t,x) that 
( |1.31c ) on [0,T] x M 3 , we can construct a local 



relativistic Maxwellian A4(n, #,u;P) (as in (1.17)) as long as the maps (1.29a) 
(1.29b) are invertible in a neighborhood of (r),p)([0, T] x M 3 ). Now for any solution 



of rE system, the identity ( 1.35 1 implies that ( 1.21 1 holds for the Boltzmann energy- 

as noted in Section [E4} the 
0. At this point 



momentum tensor constructed out of M. . Furthermore 
right-hand side of the rB equation vanishes for Ai : C(JA, M.) 
one may be tempted to conclude that M. itself is a solution to the rB equation. 
In fact, this is not the case in general; this is exactly the difficulty that Theorem 
[2] addresses. Its content can roughly be expressed as follows: if e is small enough 
and positive, then M is "close to" a solution F e of the rB equation. We remark 



that even though M. itself is not in general a solution of ( 1.1 1, it nevertheless plays 
a prominent role in our construction of a local solution of (1.1). In particular, 



the validity of the Hilbert expansion of Section |3.7| relies heavily on the fact that 
A4(n(t,x),9(t,x),u(t,x); P) is a Maxwellian and that (n(t, x), 6(t, x), u(t, x)) is a 
solution to the rE system. 



1.7. Statement of main results. In this section, we state our two independent 
theorems. The first concerns the existence of local solutions to the rE system 
satisfying the technical assumptions that are used in the proof of our main result. 
The second is our main result, which provides criteria that ensure that solutions 
to the rE system are the hydrodynamic limit of solutions to the rB equation. For 
clarity, we first state our assumptions on the initial data in a separate paragraph. 

Initial Data for the rE System. Let V = (n, 9, u 1 , u 2 , u 3 ) be initial data for 
the relativistic Euler equations under the relations (1.27a) - (1.27c), and define 



" d ~ T 2 ^"- Let f2„, z = (h, z)(]R 3 ) C (0, 00) x (0, 00) be the image of the initial data 
(h,z), and assume that there exists a compact convex set £l' nz C (0, 00) x (0, 00) 

containing Q n ,z in its interior. Let £l'„ p = (Sj, s ^)(Q.' n z ) C (0, 00) x (0, 00), where 
the functions Sj(n, z) and^{n 1 z) are defined in (1.29a) and (1.29b). Let O'' be a 



def 



compact convex subset of (0, 00) x (0, 00) containing p - Note in particular that 
our assumptions imply that there exist constants n\ > and 9\ > such that 



ni < inf n(x), 

seR 3 

9 l < inf 9(x). 
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Assume that the map (n, z) — > (Sj(n, z), CP(n, z)) == {r],p) has a well-defined, smooth 
inverse defined on Of' p that maps 0" into a compact subset of (0, oo) x (0, oo) con- 
taining Q' n z . Let p = moc 2 Ti ^|^j + 3p be as defined in ( |1.27b| , and view p, r\, and 
p as functions of (n, z) on £l' n z . Assume that on the set { (77(71, z), p(n, z)) | (n, z) € 
^nz}' P can ^ e written as a smooth function p — fkineticiV) P) an d that < 

dp 

Remark 1.2. As discussed at the end of the introduction, the careful assumptions 
on the fluid initial data are designed to ensure that they fall within the regime 
of hyperbolicity, and to ensure the invertibility of the maps between the solution 
variables (n, 9) and (77, p). However, if our two conjectures are correct, then both of 
these conditions are autom ati cally verified whenever the fluid variables are positive. 



In any case, Lemmas 3.5 and 3.6 together show that whenever 9 is uniformly small 



and positive, or in whenever 9 is uniformly large, these conditions hold. 

Remark 1.3. The additional convexity assumptions on Q! n z and Q!' p are technical 
conditions that are used in our proof of Theorem [TJ e.g. to conclude ( |2.1[ ); see 41 



Proposition B.0.4] for a discussion on the role of convexity in this context (roughly 
speaking, convexity is needed so that one can apply the mean value theorem). 

Theorem 1. (Local Existence for the rE System) Consider initial data that 
are subject to the restrictions described above. Assume that N > 3 and that there 
exist constants n > 0, 8 > such that (n, "° c 9 ) £ Q' n z (defined above), and such 
that 

(1.36) ||V-V|U r <oo, 

where V = (n, 9,0, 0, 0). Then these data launch a unique classical solution V = 
(n, 9, u , u 2 , u 3 ) to the rE system existing on a nontrivial slab [0, T] x R 3 upon which 

< inf nit, x), 

(t,at)e[0,T]xR3 

0< inf 9(t,x). 

(t,S)e[0,T]xK 3 

V has the following regularity: V-Ve C jv ~ 2 ([0, T] x R 3 ) f]^~^ C k ([0, T],H N ~ k ). 



If in addition the initial data are such that the inequalities (1.39) below are 
strictly verified att — 0, then there exists a spacetime slab [0, T"] x R 3 of existence, 



with < T' < T, upon which the condition (1.39) remains verified. 

Finally, if\\V — V||gw < 5, and S is sufficiently small, then there exist constants 
9* > and C > 0, and a slab of existence [0, T"] x R 3 , with T" > C" /S, upon which 
the following bounds are satisfied: 

(1.37) 9*<9(t,x) <26*, 

(1.38) |c-V(M)| <CS, (j = 1,2,3). 



As is shown below in Lemma \l.l\ if 5 is sufficiently small, then the bounds (1.37) 



(1.38) will also imply that the technical conditions (1.39) are verified on [0, T"] X 



3 



Remark 1.4. The conclusions of this theorem regarding the slab [0, T'] x R 3 follow 
easily from the regularity properties of the solution. 
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Remark 1.5. There are two main obstacles to extending our principal result, which 
is Theorem [2] below, to a global-in-time existence result for the rB solution. The 
first is that local solutions to the rE system tend to form shocks in finite time. In 
fact, in [IT], Christodoulou showed that there are data arbitrarily close to that of a 
uniform, quiet fluid state (i.e., V = V) that launch solutions which form shocks in 
finite time. Since our construction of a local solution to the rB equation relies on 
the availability of the solution to the rE system, the breakdown of the fluid solution 
could in principle allow for a breakdown in the Boltzmann solution. 



The second obstacle is the possible breakdown of the technical condition (1.39) 



satisfied by the fluid solution; this breakdown is sometimes avoidable. More specifi- 



cally, the condition ( 1.39 1, which plays a key role in the analysis of Section 3.8 may 



break down in finite time even before the shock happens. However, we are aware of 



a class of data that launch solutions for which (1.39) holds until the time of shock 
formation. The details are contained in we offer a quick summary. One con- 
siders initial data for the rE system that satisfy the assumptions of Theorem [T] and 
the following additional assumptions: the data are irrotational and isentropic (i.e, 
rj = const), and they coincide with the constant state V = (n, 6, 0,0,0) outside of 
the unit sphere centered at the origin of the Cauchy hypersurface {t = 0}. If the de- 
partur^J from constant state is < 5, where S is sufficiently small, then the estimates 



of 11 Theorem 13.1] imply the following fact: on the exterior of an outgoing sound 
cone C emanating from a sphere of radius 1 — e contained in {t — 0}, where e = e(<5) 
is a sufficiently small positive number satisfying < e < 1/2, n and 9 can be con- 
tinuously extendecQ to [0, T max ] x ^ 3 /Cint, where Cj„ t denotes the interior of C, and 
[0, T max ) is the maximal time interval of classical existence, i.e., the time of first 
shock formation. Furthermore, the estimates \n(t,x) — n\ < CS, \9(t,x) — 8\ < CS, 
and | C -V(i,:r)| < CS, (j = 1,2,3), hold on the region on [0,T max ] x R 3 /C mt . 
Christodoulou's theorem does not prove that the same estimates hold in Ci n t, but 
on pg. 6, he remarks that the estimates do hold on the interior, and are in fact easier 
to prove than the exterior estimates. Under these conditions, we may piece together 
the conclusions from the two regions [0,T max ] x R 3 /Ci n t and Ci n t to conclude the 



following: for sufficiently small S, an inequality of the form (1.39) is satisfied on any 
slab [0, T] x M 3 , with T < T max . Consequently, the conclusions of Theorem [2] hold 
on such a slab. 

Now that we have a large class of suitable solutions to the rE system available, we 
are ready to state our main theorem. Note that our main theorem is independent 
of Theorem [l] the role of Theorem [T] is to ensure that there are fluid solutions that 
can be used in the hypotheses of Theorem [2j 



1. 



Theorem 2. Let (n{t,x),9(t,x),u(t,xf) be a sufficiently regular (see Remark 
below) solution to the relativistic Euler equations (1.22) for (t,x) € [0,T] x ffi^ 



Construct the local Maxwellian A4(n(t, x), 9(t, x), u{t, x); P) as in \1.17). A 



that there exist constants C > 0, 9m > 0, and a € (1/2,1) such that for every 
(t,x,P) € [0,T] x El x R 3 p , the global Maxwellian J{P) = e - cp0 ^ kBeM *> from 



The notion of smallness is measured by a Sobolev norm of suitably high order. 
^Even though the L°° norm of the solution remains bounded, the Sobolev norm of the solution 
blows Up as t \ Tmax- 
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(1.16) verifies the inequalities 
J(P) 



(1.39) 



C 



< M(n(t, x), 9(t, x),u(t, x); P) < CJ a (P). 



Define initially 

6 

F e (0, x, P) = M(0, x,P) + J2 z n Fn(0, x, P) + e 3 F R , £ (0, x, P) > 0. 

n=l 

Then 3eq > such that for each < e < Eq there exists a unique classical solution 



F £ of the relativistic Boltzmann equation (1.1) of the form (1.2) for all (t,x,P) £ 
[0, T] xR| x R 3 p. Furthermore, there exists a constant Ct — Ct{M.,F\, . . . , F e ) > 
such that for all e £ (0, eq) and for any i > 9, the following estimates hold: 



e 3/2 sup 

0<t<T 



F R;s (t)/VMW) 





+ sup 


oo,i 


0<t<T 



F R , e {t)/VM{t) 



<C T <£ 



,3/2 



f r .m/VmW) 



F R . £ (0)/^M(0) 



Recall that F R . e is the remainder from ( |1.2[ ). Moreover, we have that 
sup \\F £ (t)-M(t)\\ 2 + sup ||F £ (<)-X(<)|U <C T e, 

0<t<T 0<t<T 

where the constants Ct > are independent of e. 



Remark 1.6. Conditions which would imply (1.39) are standard in the Hilbert 
expansion literature. It generally seems to be unclear at the moment how to remove 
them in the context of classical solutions. We choose to make the assumption 



(1.39) rather than a more stringent assumption of moderate temperature variation 



in order to make it clear that (1.39) is all that is needed. The condition (1.39) is 



used ensure that the local Maxwellian A4, and the terms in the Hilbert expansion 
F n , have sufficient momentum decay; as used in Section |3.8| 

Remark 1.7. We have not specified precisely the initial conditions for the terms 
Fi(0), . . ., F 6 (0), and F R , E (0) in Theorem [2] They arc constructed via the Hilbert 
expansion in Section [3.7[ after one has the fluid initial conditions as in Theorem [T] 

Remark 1.8. In the hypotheses of Theorem[2j we have not been specific in stating 
the regularity needed from the relativistic Euler equations that is required to build 
the terms Fi, . . . , Fq in the Hilbert expansion. This is because we are not attempting 
to optimize the amount of regularity needed in the Hilbert expansion. In particular, 
including additional terms in the expansion would require a smoother solution to 
the rE system. However, for the purposes of this article, it is certainly sufficient to 
have (n(t, x),6(t,x), u(t, x)) £ C 7 ([0,T] x R§) with suitable L\ integrability. 

In the next lemma, we prove that near-constant, non- vacuum fluid states neces- 



sarily verify an inequality of the form ( 1.39 1 



Lemma 1.1. Assume that the initial data for the rE system satisfy the assumptions 
of Theorem^ including the smallness assumption in 6, and let Ai(n, 9, u\ P) be the 
local Maxwellian \1.17\ corresponding to the solution. Then if 5 is sufficiently 
small, there exist a non-trivial slab of the form [0, C /S] x K|, and constants C > 0, 
M > 0, a £ (1/2,1), such that for every (t,x,P) £ [0, C'/S] x R| x R 3 p , the 
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global Maxw ellian J{P) = e- cP ° '/^ b 0m) f rom and M(n,6,u;P) verify the 

inequalities (1.39). 



Proof. To prove that ( 1.39 ) holds, we first recall that by the conclusions of Theorem 
[l] there exist a non-trivial slab [0,C'/6] x E|, and constants n m i n > 0,n max > 
0, 9* > 0, and C > 0, such that relative to our inertial coordinate system, the 
solution satisfies the following inequalities on [0, C'/S] x R| : 

< 9* < 9{t,x) < 29*, 



|<rV(i,S)| < CS, 



= 1,2,3). 



Therefore, using the Cauchy-Schwarz inequality, it is easy to show that there exist 
dimensionless constants a G (1/2, 1), C\ > 0, C*2 > 0, depending on the dimension- 



less constants ra c 
(t,x,P) e [0,C'/6] x 

(1.40) r 



'fcs^* and 5, such that if S is sufficiently small, then for all 



l 3 p we have 



i^c^dP < 



u K P H 



-'CP 



C* 2 < a _1 Ci. 



We also observe that since u, P are both future-directed and timelike, we have that 
u K P K < -\u K u K \ 1/2 \P\P A | 1/2 = -m c 2 . In particular, u K P K < 0. It thus follows 



from definitions ( |1.16 1 and (1.17), the uniform bounds on n(t,x) and 9(t,x) (which 
are strict for 9), and inequality (1.40), that for any positive constant 9m satisfying 



C 2 <m c 3 /(k B 9 M )<a- 1 C ll 
there exists a dimensionless constant C > 1, depending on the constants 



m 3 c 
holds. 



fcs^*, and 5, such that inequality (1.39) 



□ 



1.8. Historical background. There have been numerous important contributions 
to the subject of fluid dynamic limits of the Newtonian Boltzmann equation. Due 
to length constraints, it is impossible to give a comprehensive list. We only point 
out a brief few works, including the early work of Grad 24 25 . In the context of 
DiPerna-Lions 

for a more 



renormalized weak solutions, we mention the fluid limits shown 
We refer to the review articles 



33 37 



21 36 53 



in [T^|5j[T3}[22}[23 

comprehensive list of references. 

For fluid limits in the context of strong solutions to the non-relativistic Euler and 
Boltzmann equations, we mention the work of Nishida 38 , Caflisch [8], Bardos-Ukai 
j6j, Guo [2"8p9"] , Liu- Yang- Yu 34 , Guo- Jang- Jiang 31 and recently Guo-Jang [26] , 
In this work we will use the Hilbert expansion approach from Caflisch [8] combined 
with the recent developments in Guo- Jang- Jiang 31 to allow a non-zero initial 
condition for the remainder: Fn- £ (0,x, P). 

Other works which connect to and motivate different elements of our estimates/ 
results include: 11 27 43 45 47 49]- 5l]. They will be discussed in more detail at 



the appropriate time in the following developments. 

Much less is known for the relativistic Boltzmann equation. Formal fluid limit 
calculations are shown in the textbooks [9 12 . Linearized hydrodynamics are also 



studied in [15] . Our new contribution is to prove the existence of "normal solutions" 
via a Hilbert expansion to the full non-linear relativistic Boltzmann equation. It 
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would be useful to check if these types of solutions can also be constructed for the 



relativistic Landau equation, as in for instance 49 



1.9. Outline of the structure of the article. We now briefly outline the re- 
mainder of this article. In Section [2j we sketch a proof of local existence for the rE 
system that is based on well-known techniques. In Section [3j we provide additional 
background for the relativistic Boltzmann equation. In Section |3.1[ we provide an 
expression for the Boltzmann collision operator that will be used in the subsequent 
analysis. In Section [372] we study macroscopic quantities associated to a particle 
density function F(t, x, P), and we discuss the corresponding conservation laws that 



hold whenever F is a solution to the rB equation. In Section 3.3 we discuss the 
thermodynamic relations that hold between the macroscopic quantities associated 
to the special case F — A4, where A4 is a Maxwellian. In Section [3~4} we discuss 
the issue of solving for (n, 6) in terms of (77, p), assuming that the aforementioned 
thermodynamic relations hold. In Section |3.5| we prove that there are regimes of 
hyperbolicity for the rE system under the same thermodynamic relations. In Sec- 
tion |3.7[ we carry out the Hilbert expansion for F e in detail. Finally, in Section 
|3.8| we provide a proof of our main main theorem. 



2. Local existence for the relativistic Euler equations 

In this section we sketch a proof of local existence for the rE system; i.e, we 
sketch a proof of Theorem [T] of Section |1.7| During the course of the proof, which 
is mostly along standard lines, we provide references that indicate where one may 
find the omitted details. 



Sketch of the proof of Theorem [TJ Theorem [T] can be proved using the 
method of energy currents, a technique which was first applied to the rE system 
by Christodoulou in 



11 



43 , which 



For complete details, one can first look in 
describes in detail how to derive energy estimates for the linearized rE system; we 
sketch this derivation below. Once one has these energy estimates, the proof of 



local existence for the non-relativistic Euler equations given in 35 can be easily 

for the 



modified to handle the rE system. One may also consult 32 , 40 



or 



43 



essential ideas on how to finish the proof once one has energy estimates for the 
linearized system. Furthermore, although we do not need it for this article, we 



remark that 43 contains a proof of continuous dependence on initial data. 

In order to apply this method, we use , p, u , u 2 , u 3 ) as our unknowns for 

the equations ( 1.31a I - (1.31e). The assumptions on the data allow us to use the 



fun ctions fj and *}3 from (|1.29a 

1 „°,2 



1.31 V d = (h, 



(1.29b) to transform the initial data (see Remark 
= f (f],p, u 1 , u 2 , u 3 ) such that 



to initial data W 



(2.1) 

and such that 



|W- W|| H w < oo, 



< inf fi(x), 

< inf p(x), 

sm 3 
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where W = (77, p, 0, 0, 0). Here, r\ and p are equal to $)(n, 0) and ty(ri, 6) respectively. 



e.g. the appendix of 43 



That (2.1) follows from (1.36) can be shown via Sobolev-Moser type estimates; see 



Typical proofs of local existence are based on cither the construction of con- 
vergent sequence of iterates, or a contraction mapping argument. Both of these 
arguments require that one prove energy estimates for the linearized Euler equa- 
tions, which are the following system: 

(2.2a) uPdrf = & 



(2.2b) 
(2.2c) 



u^d^p + q^dou* + qd k u k = ©, 



Mk . 

~ def / 1 ~2 ~3\ 



Cj = 1,2,3). 



In the above equations, W = (i],p, u 1 , u 2 , u 3 ) is the "background" (which can be 
thought of as the previous iterate in an iteration scheme), W = f (f],p, ii 1 , u 2 , ii 3 ) 
is the variation (which can be thought of as either the next iterate or one of its 

spatial derivatives), and the terms b = (^,®,Sj 1 ,Sj 2 ,^) 3 ) are the inhomogeneous 
terms that arise from the iteration + differentiation procedure. Furthermore, u° = 

\J c2 + Tfj=i(^ J ) 2 ; = u»u v + (g~ l y u , and q is the function of the state-space 



variables from (1.31e) evaluated at the background. 



To deduce energy estimates for the linearized systems, one first defines energy 
currents 3, which are vectorfields that depend quadratically on the variations W: 

(u k u k ) 2 - 



(2.3) 3°(W,W) 



u° . 2 n u k ii k . 
rf + —p 2 + 2-^-p+(p+p)u L 
q u v 



3° (W, W) = v7i] 2 + ^rp 2 + 2u>p + (p + pju 1 



• k ■ 

u u k 



u k Uk 



{u k u k ) 2 



(u ) 2 



U = 1,2,3). 



The two key properties of the above energy currents, which are discussed in detail 
are (i) it can be shown that <j°(W, W) is a positive definite quadratic form 



43 



in W, and (ii) if W is a solution to the linearized rE system (2.2a) - (2.2c), then it 
can be shown that 9 K [3 K (W, W)] does not depend on the derivatives of W. More 
specifically, the following formula holds: 



(2.4) M 0" = (0 /t «'V + / ,ly ]1>~ 



ii k u k 



(UkU 



k\2 



2u k u k {p + p) 



2f)$ + 2- 



where 6, -f) 1 , -f) 2 , -f) 3 ) are the inhomogeneous terms in (2.2a|- (2.2c) 
We therefore can define an energy £(rj) > by 



(2.5) 



£ 2 (i) 



a°(W,W) da, 



and £ can be used to control ||W||/,2. More specifically, property (i) can be used to 
show that if q, p, p are uniformly bounded from above and away from on [0, t] x IR 3 , 
and the are bounded from above on [0,t] x M 3 , then there exists a constant 
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C > depending only on the values of W, || c?t W|| , and ||9xW||i«., such that 
on [0,t], we have 



(2.6) 



C- 1 ||W|| i2 < £ < C||W|| i2 . 



Remark 2.1. If q < 0, then 3° is no longer positive definite, and inequality (2.6) 
fails. Since q 



„2 dp 

- dp 



(p + p), it is clear that the non- negativity of 



Op 



plays a 



fundamental role in the well-posedness of the rE system. This explains the signif- 
icance of Lemma |3.6| and Conjecture [2] Furthermore, we note that the conditions 
on the initial data guarantee that q[x) = q(0,x) is uniformly positive. 



Furthermore, the Cauchy-Schwarz inequality for integrals, (2.4), and (2.6) imply 
that 



(2.7) 



||W(W,W)]|| L1 < C||W|| L 2||b|| L 2 < C£||b|| L ». 



Using (2.5), the divergence theorem (assuming suitable fall-off conditions at infin- 
ity), and (2.7), we have that 



d K [3 K (W,W)] dx 
<C£||b|| za , 



which leads to the following energy estimate for the linearized system: 



(2.8) 



-£(t) < C||b|| i2 . 



The availability of (2. 



is the fundamental reason that the rE system has local- 
in-time solutions belonging to a Sobolev space. This concludes our abbreviated 
discussion of the existence aspect of Theorem [T] in terms of W; we comment on the 
variables V = (n, 0, u l , u 2 , u 3 ) near the end of the proof. 

Let us now assume that we have a local solution W to the nonlinear rE equations 
near the constant state W; we will make a few remarks about the time of existence. 
In the nonlinear case, one can define energies (with N > 3 so that various Sobolev- 
Moser type inequalities are valid) 



£ 2 



N 



{t) d = f 



E 

0<|5|<Af 



/ 3°{d s (W- W),3 S (W- W))dx. 



Furthermore, all of the inhomogeneities that arise upon differentiating the rE equa- 
tions are of quadratic order or higher. Consequently, near a constant state, the 



inhomogeneous terms analogous to the term ||b||^2 on the right-hand side of (2.8) 
can be bounded by C£, 2 N {t), and the resulting energy estimate is 

d 
Jt 



(2.9) 



^M*) < ce 2 N (t). 



Consequently, we may apply Gronwall's inequality to (2.9) and use property (2.6) 
to deduce an a-priori estimate of the form 
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It follows from (2.10) that the time of existence is at least of size C/S. This fact is 
a consequence of a standard continuation principle that is available for hyperbolic 



PDEs (consult e.g. 32 Chapter 6] or [42] for the essential ideas) , which implies that 
the solution exists as long as the a-priori energy estimates lead to the conclusion 
that II W — W||jjjv is sufficiently small. 



We remark that in the formula (2.10), C is a numerical constant that depends 
on an a-priori assumption concerning the subset of (0, oo) x (0, oo) to which the 
pair {rj,p) belongs, and on B, where B > is a fixed a-priori upper bound for 
|W — W||jjw. These a-priori assumptions imply that the solution escapes neither 



the regime of hyperbolicity, nor a convex subset (see Remark 1.3) of the regime in 
which the map (n, z) — » \S)(n, z),ty(n, z)j = (r],p) is invertible; on the time interval 
[0, C/S\, these a-priori assumptions can be shown to hold through a bootstrap ar- 
gument. Furthermore, once we have shown (2.10), this inequality can be translated 



into an inequality for the original variables V; Sobolev-Moser type estimates allow 
us to estimate (here we again use the convexity assumption discussed in Remark 



1.3) ||V — V|| ff « < C||W — W|| H jy. This is possible because V can be written as a 



smooth function of W whenever (rj, p) belongs to the domain of the inverse of the 
aforementioned map. 

3. The relativistic Boltzmann equation 

In this section, we provide some additional background material on the rela- 
tivistic Boltzmann equation. The rB equation is often expressed in the physics 



literature (see e.g. |7[ |12| ) in the Lorentz-invariant form of ( 1.1 1 where the collision 
kernel C is defined by 



C(F,G) 



dQ f dQ' f dP' 



2 ./»3 Q° 7»3 Q'° 7m P'° 



W(P, Q\P', Q')[F{P')G{Q') - F(P)G(Q)} 



In our analysis below, it will often be convenient to divide both sides of ( |1.1[ ) by 
P . We therefore introduce the normalized velocity P, defined by 



P def _c_ p 



and fi(.,.) 
(3.1) 



ttC(-, ■)■ Then (1.1) is equivalent to 



d t F 



P ■ d s F = -Q(F, F), 

£ 



where P ■ d s = P 1 d 1 + P 2 d 2 + P 3 d 3 . 

On the right-hand side of the above expression for C, the variables P and Q 
represent the pre-collisional four-momenta of a pair of particles, while P' and Q' 
represent their post-collisional four-momenta. We assume that the particle colli- 
sions are elastic, in which case the conservation of energjj^] and 3-momentum can 
be expressed as 

(3.2) P^ + Q" =P^ + Q^, (^ = 0,1,2,3). 

The transition rate W(P,Q\P' ,Q'), which is a Lorentz scalar, can be expressed as 

(3.3) W(P,Q\P',Q') = sa(g,'d)6^(P fl + Q" - - Q ,fl ), 



In our inertial coordinate system, the "energy" of a four-vector is c times its component, while 
its 3-momentum comprises its final 3 components. 



HILBERT EXPANSION FROM RELATIVISTIC BOLTZMANN TO EULER 



2:S 



where 5^ is a Dirac delta function, and a(g,'&) is the differential cross-section or 
scattering kernel. The other quantities are defined in (1.121, (1.131, and (1.14). 
We have several remarks to make. First, the Lorentz invariancc of the left-hand 



side of (1.1 1 is manifest, while the Lorentz invariance of the right-hand side follows 



from that of j£o (and similarly for Q,Q'); see (1.9). Next, we remark that i? in 



(1.13) is well-defined under (3.2), but that it may not be in general [181; i.e., in 



general, the right-hand side of ( 1.13 1 may be larger than 1 in magnitude. Finally, we 



point out that both the left and right-hand sides of (3.1) are functions of (t,x,P) 



3.1. Expression for the collision operator. In this section, we provide an alter- 



nate expression (3.4) for the collision operator, which is derived from carrying out 



certain integrations in a center-of-momentum frame. This is the expression for the 
collision operator that we will use in our analysis for the remainder of the article. 
We remark that yet another expression for the collision operator was derived in 
20| ; see 46 47 for an explanation of the connection between the expression from 



20 



and the one in (3.4) below 



One may use Lorentz transformations as described in 12 and in some detail in 



46 48 to reduce the delta functions in (3.3), thereby obtaining 



(3.4) 



Q(F,G) 



dQ / dto v a(g,d) [F{P')G{Q') - F(P)G{Q% 



where v = v (P,Q), the M0ller velocity, is given by (1.11) 



The post-collisional 3-momenta in the expression (|3.4|) can be written as follows: 
P + t 



P' 



(3.5) 



_ + |(„ +(7 _ 1)( p +0) (^ 

^_|( B + (7 _ 1)(P + 0) C^ 



where 7 = (P° + Q )/^, and • is the ordinary Euclidean dot product in M 3 . The 
energies can be expressed as 

P + Q 



P" 



Q'° = 



2^s 



=u-(P + Q), 



P° + Q° 



2^s 



u)-(P + Q). 



It is clear that P,Q, P' ,Q' satisfy (3.2). Additionally, it is explained in 48 that 



after carrying out the integrations in a center-of-momentum frame, the scattering 
angle becomes a function of P, Q, and satisfies cos?? = k ■ w, where uj g § 2 C K 3 , 
and k = k(P, Q) 6 S 2 C R 3 . We remark that an expression for k(P, Q) is given in 



Equation (5.37)] and 46 , but that its precise form is not needed here. 



3.2. Macroscopic quantities and conservation laws for rB. In this section, 
we study macroscopic quantities associated to a particle density function F(t, x, P). 
We also discuss conservation laws that hold whenever F(t, x, P) is a solution to the 
relativistic Boltzmann equation (1.1). This material is quite standard [12], but we 



include it for convenience. We begin by recalling that the energy- momentum tensor 
TjBoltz [-^1 an d the particle current Ig oltz [F] for the relativistic Boltzmann equation 
are defined in (1.20). Note that the term ^ on the right-hand side of (1.20) is 



exactly the canonical measure defined in (1.7). 
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We now define the following macroscopic quantities: 



dP 



(3.6a) S» d - f ~k B cJ s P»F(t, x, P){Hh 3 F(t, x, P)} - l} pn . 



(3.6b) -c 2 n 2 = g KX I% oltz I^ oltz , 



def m _l ; 

' Boltz' 



(3.6c) 

(3.6d) u K u K = -c 2 , 

(3.6e) p = cr 2 u K u x T^ ltz , 

(3.6f) P ^ f n^ to ir A , 

(3.6g) 77 ^' -c" 2 ^ 1 ^. 

Here, 5 is the entropy four flow, n is the proper number density, u is the four- 
velocity, p is the proper energy density, p^ v is the pressure tensor, rj is the entropy 



per particle, and IP*" c = c~ 2 u fJ, u L ' + (g^ 1 )^ is as defined in (1.26). 

We also decompose the pressure tensor by defining the pressure p and the viscous 
pressure tensor v^ v : 



p^^pU^ + v^, 



where 



(3-7) p d ^\n KX p KX = l -Xl^T£ ltz . 

We remark that H kX v kX = 0. 

The following lemma will be useful in verifying numerous identities, especially 



those of Proposition 3.3 



Lemma 3.1. Let A^ y be a matrix (of real numbers) such that 

A K ^A\g KX = g^, (p, v = 0, 1, 2, 3), 

and such that det(A) = 1 (i.e. a proper Lorentz transformation). Consider the 
change of coordinates P 4— P on T X M from one inertial coordinate system to 

def — 

another induced by A : P M = A^ K P K . Then under this change of coordinates, 1^ , u^, 
and transform as the components of a vector, while Tg U oltz and p^ v transform as 
the components of a two-contravariant tensor. In an arbitrary coordinate system, 
these statements take the following form: 



(3.8a) 


I" 


= A" W J", 


(A* = 


0,1,2,3), 


(3.8b) 


u» 


= A^ K , 


(/* = 


0,1,2,3), 


(3.8c) 


S n 


= A" K S K , 


(/* = 


0,1,2,3), 


(3.8d) 


1 Boltz 


— yv k iv \ L Boltzi 




-0,1,2,3), 


(3.8e) 


p^ 


= A» K A\r\ 




-0,1,2,3). 



Proof. We give only the proof of (3.8a); the remaining statements follow simi- 
larly. Equation (3.8a) is equivalent to the following statement, where we abbreviate 
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(3.9) 



f P»F(P l ,P 2 ,P 3 )^= [ A^P K F(A 1 a P a ,A 2 l3 P^ 1 A 3 1 P^)- 

JR 3 " JPfR 3 



F(t,x,P) = F(P) 

) dP 

I PGR 3 P° 

Recall that Lorentz transformations preserve the form of the metric g, so that 
fl^v = 9ftv = diag(— 1, 1, 1, 1), which implies in particular that P° = \J m\c 2 + \P\ 2 

and P a = \J tuqC 2 + \P\ 2 . It thus follows from the discussion in Section |l.2[ 



and 



in particular equation (1.9), that ^ 



dP 
P o 



Equation (3.9) now follows from the 



standard change of variables formula for integration. 



□ 



The conservation laws for the relativistic Boltzmann equation are given in the 
following lemma. 



Lemma 3.2. |9j Chapter 2] Let F(t,x,P) be a C 1 solution to (1.1). Then the 
conservation laws (1.21) hold. 



Note that (1.21) corresponds to the fluid conservation laws (1.22). 

3.3. Macroscopic quantities for a Maxwellian A4 = M.(n, 9, u; P). In this 
section, we prove the following proposition: 

Proposition 3.3. Let n(t,x) > 0, 9(t,x) > be positive functions on M, and let 
u(t,x) be a future- directed vectorfield on M satisfying (1.24). Consider the corre- 
sponding local relativistic Maxwellian A4(n,9,u; P) defined in (1.17) with (1.18). 
Then the following relations hold for the quantities defined in (3.6a) - (3.7): 

(3.10a) 



(3.10b) 
(3.10c) 



I»[M] = n[MK, Gu = 0,l,2,3), 
p[M] = k B n[M}9 = m c 



2 n[M] 



p[M] = 3p[M] + m c 2 n[M] 



z[M]' 
Ki(z[M]) 
K 2 (z[M}) 

p[M] 



(3-10d) TZ ltz [M]=TZJM], (/*,!/ = 0,1, 2, 3), 



(3.10e) 



\M} = 



' fluid 

4:7re A mlc 3 h~ 3 exp 



\k R J 



4irmlc 3 h 3 exp 



kr 



r,^ K 2 {z[M] 
z[M] 
r,\ K 2 (z[M }) 
z[M] 



exp 



(z[M] 



2{z[M\), 

KMM]) \ 
K 2 {z[M])J- 



By Tf luili [A4], we mean the energy-momentum tensor that results from inserting 
p[Ai] and p[M] into the expression on the right-hand side of (1.23). 

Proof. It is well-known that since u(x) is a future-directed vector satisfying ( |1.24 1, 
there exists an inertial coordinate system in which (u°, u , u 2 , u 3 ) = (c, 0,0, 0) (at 
the spacetime point x); such a frame is known as a "rest frame" for u. By Lemma 



3.1 it suffices to check that (3.10a) - (3.10c) hold in such a rest frame. Using 
~ , it foil 



definition (1.20), it follows that I ^ can be expressed in a rest frame for u as 

30 ^ dP 



(3.11) 



47ttoqC 3 K 2 (z) 



7 R3 V m c / 



per 
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By symmetry, it follows that P*[.M] is proportional to (1,0,0,0), and therefore 
also to u. Let us denote the proportionality constant by A, so that in any inertial 
coordinate system, we have that 

I»[M] = Au^. 



We thus have that A = I°/u°. Furthermore, it follows from (3.11) that in a rest 
frame for it, we have that 



(3.12) 



A 



4irm%c 3 K 2 (z) J R3 



exp 



/-zP°\ 



dP. 



We carry out the integration in (3.12) using spherical coordinates for P, and making 
the change of variables A = zP° / \jiiqc). This implies that |P| = mocz _1 (A 2 — z 2 ) 1 / 2 , 
d\P\ = to cz _1 A(A 2 — z 2 ) -1 / 2 , and we obtain 



(3.13) 



3„3„-3 



A = n- ., AiTmlc A z 

A-Kmlc 3 K 2 {z) 



Xe- x (X 2 -z 2 ) 1/2 d\ 



\=z 



From the Bessel function identity (3.44) in the case j = 2, it follows that A = n, 
which completes the proof of (3.10a). 

To prove (3.10b), we note that in a rest frame for u, it follows that P K P X H K \ = 
\P\ 2 , where is defined in (1.26). Inserting this formula into definition (3.7), 
using definition (1.20), integrating with spherical coordinates for P, and using the 

-zP°\ dP 



integration variable A as above, we have that 



(3.14) 



p[M] 



1 z 

— CTI 

3 4irm 3 c 3 K 2 (z) 
1 z 

— CTI 

3 4Trm 3 c 3 K 2 (z) 



|P| z exp 



. 4 4 —4 

47rm c z 



(— ) 



m ( 

oo 

e" 



P 
A (A 2 - 



, 2 ) 3/2 dA. 



-, which proves ( 3.10b[ ). 
To prove (3.10c), we note that in a rest frame for u, u K u\P K P x = c 2 (P ) 2 . 
Inserting this formula into definition (3.6e), carrying out the integration in spherical 
P coordinates, using the change of variables A as above, referring to definition 
using 

p[M] = : ''a - 



Referring to definition (3.43), it follows that p[M] 



(3.43), and using the recursion formula (3.45) in the case j = 2, we have that 

z 



(3.15) 



cn 



AttttIqC 3 K 2 (z) 

z 

4Trm,QC 3 K 2 (z) 

z 

4nmQC 3 K 2 (z) 



P°exp 



Amnio 3 K 2 (z 



Airnikc^z 4 



a 4 4 —4 

47rm n c z 



-Airmkc A z 



/-zP°\ 

V TTinC / 



dP 



TTJoC 
A— oo 

\=z 
A— oo 

A=z 



e - A A 2 (A 2 -z 2 ) 1 / 2 dA 
e - A (A 2 -z 2 ) 3 / 2 dA 
e- A (A 2 -z 2 ) 1 / 2 dA 



TOq c n 



■3K 2 (z) , K^z) 



K 2 (z) 



SnksO + m c 2 n }\ \ - 2 

K 2 (z) 



"^oc 2 n[7W]|^4 " k B n[M}9. 
K 2 {z) 



Combining this identity with (3.10b), we deduce (3.10c). 
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To prove (3.10d), we first note that in a rest frame for u, we have that 
T&JM] = di ag (p[M},p[M},p[M],p[M]), 

and 



(3.16) 



TZtAM]=cn 



Awmlc 3 K 2 (z 



[ P^P V expf 



-zP°\dP 
itiqc J P° 



The fact that /j / 1/ =>■ T£\ t Z [M] — follows from symmetry. The fact that 
T B ° oltz [M] = p[M] follows directly from comparing the integral expressions (3.15) 
and ( |3 . 16 1, while the fact that T B J oltz [M] = p[M] (there is no summation in j here) 
follows from comparing the integral expressions ( 3.14[ ) and (3.16) using symmetry. 

To prove (3.10e), we notice that in a rest frame for u, we have that u K P K = —cP°. 
Inserting this formula into definitions (3.6a) and (3.6g), and evaluating the two 
integrals that arise as in (3.13) and (3.15), we have that 

-zP^ 



r)[M] = {-k B i 



-zP° 

TUqC 



47ttoqC 3 K 2 (z) 



In 



exp 



(— ) 



= k f 



z 



exp 



-zP° 



fcB 47rm3 C 3^ 2 ( 2 )l ln 

z 

47rmQC 3 K 2 (z) tuqC 



m c 

z 



m c 
l\dP 

P°dP 



AttitiqC 3 K 2 (z) 



]-'L 




f =p( _iP °) 




R3 V m Q c J 



dp 



3K 2 (z) , Ki{z) 



kB 4Trm 3 c 3 K 2 {z){ 
(3.17) =kB {3 + z X&}- kB { 



h 3 n 



z 



In 



■iwrriQC 3 K 2 (z) 

,? z 
h 3 n 



1 \4irml c 3 



K 2 {z) 



AttitIqC 3 K 2 (z) 



(3.10e) now follows from ( |3.17 l and simple algebraic manipulation (solve for n). □ 

The next lemma was used in Section 2] to derive an equivalent version of the 
rE system; i.e., equations (1.31a) - (1.31e). More specifically, only equation (3.19) 
was used. However, as an aside, we also discuss the fundamental thermodynamic 
relation (3.18) (see 10 ), which can, in consideration of the positivity of n and 9, 
be used to show that rj can be written as a smooth function of n, p. 



Proposition 3.4. Assume that the functional relations (3.10b), (3.10c), and (3.10c) 
hold for the variables p, p, n, ( 



(3.18) 
(3.19) 



and rj. Then the additional relations also hold: 
dp 



n9 



dr) 



p + p = n 



dp 
On 



Proof. To ease the notation, we use the notation ( 1.18 ) and abbreviate Kj — Kj (z), 
K'j = -^Kj(z), where Kj(z) is the Bessel function defined in (3.43). To begin the 
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proof of (3.18), we first note that by the chain rule, it follows that 
(3.20) 



dp 


dp 


dz 


drj 


n ~ dz 


n 9 V 



We claim that equation (|3.20|) leads to the following identities: 



(3.21) 



dp 
drj 



nrriQC 
k B z 



A' 



K^K' 2 K 3 



+ z- 



z~ 2 K 2 + K' 



K^K' 2 K 3 • 



n6. 



which completes the proof of (3.18). 



To see that (3.21 ) holds, we note that differentiating the last equality in (3.10c) 



leads to the relation 
(3.22) 



dp 

dz 



K 2 K 3 



while differentiating each side of (3.10e) with respect to r\ (while n is held constant) 
leads to the following identity: 

dz 



(3.23) 



dr] 



K' 3 



l K 3 



K^K' 2 K 3 



k B z 



- X K 2 {z). 



Inserting (3.22) and (3.23) into the left-hand side of (3.21 ) implies the first equality. 
The second equality in (3.21) follows from (3.46) in the case j — 2, which implies 
that the term above the under-braces is equal to 1, and from the definition of z. 



The proof of (3.19) follows similarly using the chain rule identity 



dp 

c) 1 1 



_i_ dp 
^ dz 



1 , and we omit the calculations. 

cm I n ' 



□ 



3.4. The invertibility of the maps Sj(n,z) and Cp(n,z). In this short section, 
we state and prove Lemma 3.5 which addresses the issue of solving for (n,0) in 
terms of (77, p). This lemma rigorously shows that Conjecture [l] is true outside of 
a compact set of 9 values. Furthermore, at the end of the section, we provide 
Figure [T] which is our numerical evidence for the validity of the conjecture. For the 
purposes of avoiding repetition, during the proof of Lemma |3.5[ it is convenient to 
use notation that is defined below in the proof of Lemma |3.6| However, logically 
speaking, the proof of Lemma |3. 5 1 comes before the proof of Lemma |3.6| 



Lemma 3.5. Consider the smooth maps 77 = Sj(n,z) and p = CP(n, z) defined 
(1.29a) and (1.29b) respectively. Then the map (n,z) — > (f)(n, z), CP(n, z)) is in- 
vertible with smooth inverse if0<z< 1/10 or z > 70. 

Proof. Using (1.27a) and (1.27c), it follows that pexp (^J is a smooth function 
of z alone. Therefore, by the implicit function theorem, since p > holds whenever 
z > 0, we can locally solve for z in terms of 7/, p if ^ 

< holds for < z < 1/10 and z > 70. 



7^ 0. We will show that 



()p 



(3.24) 



We begin by quoting equation (3.31), which states that 
d z \ v p 



K 2 (z) \K 2 (z) 



Ki(z) 
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Using (3.24) and the expansions (3.331 - (3.34), it follows that 



Thus, using the bounds (3.35), we obtain that 

*M£ < -3, 



(0 < z < 1/10). 



Since p > 0, it follows that d z \„p whenever < z < 1/10 as desired. 

On the other hand, using the expansions ( 3.38| ) and (3.39), it follows that 



z d AR = _ 5/2 + ^ + 3z7 x (z) + z%(z). 
p 8z 



Therefore, using the bounds (3.40), it follows that 
dz\ v P 



P 



< -1, 



(z > 70). 



Therefore, by the implicit function theorem, we can solve for z in terms of 77, p if 
< z < 1/10 or z > 70. Since n — m,Q 1 c~ 2 pz from (1.27a), the same is true of n. 

□ 



This completes the proof of Lemma [3. 5 1 
Remark 3.1. It is clear from the proof of the lemma that Conjecture [I] can be 



Op 



for all z > 0. Thus, the numerical 



shown by demonstrating the negativity of ^ 

plot in Figure [T] which was created with Maple 11.0, is the motivation for our 
conjecture. 

3.5. Regimes of hyperbolicity for the rE system and the existence of the 
kinetic equation of state. In this section, we prove that whenever 9 is sufficiently 
small and positive or sufficiently large, there exists an equation of state of the form 
(1.28), i.e., of the formp = fkineticiVi p)- Furthermore, under the same temperature 
assumptions, we show that the equation of state satisfies < a ^ kinc 



dp 



< 1. We 



remark that 



dfk 



dp 



can be expressed as a function of 9 alone. As previously 



discussed, this condition is sufficient to ensure the hyperbolicity of the rE system in 
these temperature regimes; in particular, as discussed in Remark |2.1[ the condition 



< 



9fki- 



Dp 



plays a fundamental role in the proof of local existence. Our result 

rigorously shows that outside of a compact set of 9 values, a slightly weaker version 
of Conjecture [2] holds. Furthermore, at the end of this section, we provide Figure 
[2] which is our numerical evidence for the validity of the conjecture. Also see the 
discussion at the end of Section |1.6| For convenience, we use the variable z from 
( |1.18 1 during the statement and proof of the lemma. 

Lemma 3.6. (Hyperbolicity of the rE system) Assume that the functional 
relations (3.10b), (3.10c), and (3.10c) hold for the macroscopic variables n,9,r],p, 



and p. Then if < z < 1/10 or z > 70, p can be expressed as a smooth function 

f kinetic Of 77 and p : p = fkineUciW, P)- 

Furthermore, the following estimate holds for < z < 1/10 : 

dp I 1 

by-n ~ 3 



(3.25) 



< z z 



30 



JARED SPECK AND ROBERT M. STRAIN 



Figure 1. z b 



plotted as a function of z. 



Additionally, the following estimate holds for z > 70 : 



(3.26) 



dp , 3z 
dp^~ 5 



< 41. 



Remark 3.2. We did not attempt to be optimal in our estimate of the error terms 
on the right-hand sides of the above inequalities. 



Proof. It follows from (|3.10b|), (|3.10c|), and (|3.10e|) that 
(3.27) 



p = 4ne 4 mQC 5 h 3 exp 



-V\K 2 (z) ( Kx(zY 



( :5.2m /, = / '( z §f+ 3 



We use the following version of the chain rule: 
further deduce that 



%\ v = Sf- Usin s S we 



(3.29) 



d z \ v P 



= 3 + z- 



-?-Y 



K 2 (z) \d z \ vP Jdzl K 2 (z) 
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Using the bounds (3.40) and the expression (3.41 ), it can be checked that for z > 70, 
we have that 



(3.42) 



d z 



nP 



0, 



3z 



< 41. 



The facts that p can be expressed as a smooth function of 77 and p whenever z > 70, 
and that inequality (3.26) is verified, both easily follow from (3.42). □ 



Remark 3.3. Notice that the Conjecture [2] is equivalent to the conjecture that the 
right-hand side of (3.32) is > 3 for all z > 0. In Figure [2] we present a numerical 



plot, which was created with Maple 11.0, that covers the set of z values lying outside 
of the scope of Lemma [3~6l and that suggests that this conjecture is true. Note that 
the inequalities of the conjecture are stronger than those proved in the lemma. 




Figure 2. 



plotted as a function of z. 



3.6. Bessel function identities and inequalities. We now state the technical 
lemma that contains the Bessel function properties that we have used throughout 



this article. The expansion (3.47) (including the error terms) and inequality (3.48) 
can be found in 39 . The remaining identities can be found in [12] Chapter 2]. 
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Lemma 3.7. (Properties of Bessel functions) Let Kj(z) be the Bessel function 
defined by 



Then the following identities hold: 



- A (A 2 - z y~w d\, (j>o). 



(3.44) K 3 {z) = 



1)! 1 

(2j-2)i ~Z j x _ 



A— oo 



Ae- A (A 2 -^ 2 )^( 3 / 2 )dA, (j>0), 



(3.45) X J+1 (z) = 2j^^ + i^-i(z), (j > 1), 



also 



(3.47) 



(i > o), 

/ Tl— 1 

" m=0 



{] > 0, n > 1) 



where the following additional identities and inequalities also hold 
A 3 ,o = 1, 

(4j 2 -l)(4j 2 -3 2 )---(4j 2 -(2m-l) 2 ) 



A, 



| 7 ,,„(z)| < 2exp([j 2 - l/^z- 1 )!^!, (j > 0,n > 1), 
(3.48) /^-(z) <i^ +1 (z), (j>0). 



, (j > 0,771 > 1) 



The following corollary of Lemma |3.7| is used in the proof of Lemma 3.6 
Corollary 3.8. For < z < 1/10, the following inequalities hold: 

Kx(z) z 



K 2 {z) 2 
Ki{z)\2 z 2 



\K 2 (z)J 



(3.49a) 
(3.49b) 

For z > 10, the following inequalities hold: 
(3.50a) 

(3.50b) 



< 2z 2 , 



< 2z J 



Ki(z) 3 15 

K 2 (z) 2z 8z 2 



z z z 



\K*{z)) 



16 

Z J 

40 
< — . 

- r3 



Proof. We remark that throughout the proof, we make no attempt to be optimal in 
our estimates. Using (3.43) in the case j = 1, and the fact that f\~£° Ae~ A dX = 1, 
it follows that 



zKx{z) - 1 



e A dX 



A=0 



Ae" 



1 



1 \dX. 



The first integral is trivially bounded in magnitude by z 2 /2. Using the fact that 
I a/ 1 — (z/A) 2 — l| < (z/A) 2 on the domain < z/A < 1, it follows that the second 
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integral is bounded in magnitude by 



A=0 



e~ x d\ < z. 



We therefore conclude that 
(3.51) 



K x (z)- 



< 1 



Using (3.44| in the case j — 2 and similar arguments, which we leave to the reader, 
we also conclude that 

(3.52) 



z 

^ 1+ 3- 



Using (3.51) and (3.52 1, together with simple algebraic estimates, it follows that 



for < z < 1/10, we have 



K 2 (z) 



< 2z* 



This proves (3.49a). Inequality (3.49b) follows from similar reasoning; we leave the 
details to the reader. 

To prove ( 3.50a[ ), we first decompose 

Kx{z) _1 + A 



(3.53) 

where 
(3.54a) 

(3.54b) 



l + B 

A = 
B = 



(1- 




- B + B 2 


3 


15 


, 71,3 


8^ " 


128z 2 


3 ' 
Z 6 


15 


105 


72,3 


8^ " 


V 128z 2 


' 3 ' 



B 3 
l + B 



Y 



and the 7j jn are from (3.47). 

For the remainder of the proof, we will now assume that z > 10; all of our 
estimates will hold on this domain. Now using (3.47), it can be checked that the 
following inequalities hold: 

(3-55) | 7 i, 3 | < \, |72, 3 | < 1. 



Consequently, it is easy to check that the following estimates hold: 
(3.56a) 



(3.56b) 
(3.56c) 
(3.56d) 



1 , 1 

— < A < — , 

4z ~ ~ 2z 

1 2 

-<B<-, 
z z 



B 



15 



< 



1 



B — [ — 
82 



15\ 2 



4 

< — . 

- ^3 



Using simple algebraic calculations, it follows from the expansions (3.53), (3.54a), 
and ( |3.54b| > that 

K x {z) =1 _ 3_ ^5 
2z 8z 2 



(3.57) 



K 2 (z) 



0{z~ 
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In (3.571, the symbol 0{z~ 3 ) denotes the cubic (in z~ l ) and higher-order terms 
that arise in the expansion of W . We now estimate this 0(z~ 3 ) term by using 
the expansions (3.531 - (3.54b) to split it into the following 3 pieces: 

B 3 



(3.58a) 
(3.58b) 
(3.58c) 



/= - 



15 

128z 2 



II 
III 



B 



:{"( 



71,3 72,3 



71,3 
Z 3 

15 

82 



B 2 



- B + B z - 



l + B 



B 2 



B 3 



l + B 



}■ 



15\ 2 

8zy 



B 3 



l + B 



It is easy to see by sign considerations (i.e., using B > 0) that | — B + B 2 — jqrg | < 
\B\. Using also (3.55) and (3.56b), we conclude that the following inequality holds: 



(3.59) 



l/l < 



15 



71,3 



128z 2 



\B\ < 



15 



128z : 



4z 2 



2 3 

- < . 

z 4z 3 



For the term II, we use similar sign considerations, together with the estimates 



(3.56b) and (3.56c) to conclude that 
3 



(3.60) 



\II\ < 



8z 



B 



15 



Finally, for the term III, we use the fact that B > 0, together with (3.55), (3.56b I, 
and ( |3.56d 1 to conclude that 



(3.61) 
\HI\ < { 



7l,3 


+ 


72,3 


+ 


B 2 — 


Q 2 


z 3 




z 3 







\ b \ 3 }< A+- 

1 1 J - 4z 3 z c 



4 8 53 
I-' 



Adding (3.59), (3.60), and (3.61 1, we arrive at (3.50a|. 



Inequality (3.50b) can be shown directly from (3.50a); we omit the details. □ 



3.7. The Hilbert expansion. In this section, we perform a Hilbert expansion 
for the rB equation (3.1). We decompose the solution F £ as the sum (1.2) where 



F , Fi, . . . , Fq in (1.2) will be independent of e. Also, -Fr ;e is called the remainder 
term; it will depend upon e. Our main goal in this section is to explain how one 
can prove Proposition |3.9[ which summarizes the behavior of Fq, F\, . . . , Fq\ the 
remainder term F^ £ is analyzed in detail in the next section. 
We begin by inserting the expansion (1.2) into (3.1) to obtain 



k=0 



J2e k (d t + P-d s )F k +e 3 (d t + P-d s )F R:£ = J2 E ^^Q^Fj) 

fc=0 i+j = k 
0<i,j<6 

6 

+ s 5 Q(F R;e ,F R , e ) + J2 £ 2+k {Q(Fr; S , F k ) + Q(F k , F R , £ )} + A. 



fe=0 
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Above A = f 



Efcly J2i+j=k,i<i,j<6 £l+J 1 Q{ F ^ F j)- Equating like powers of e on 
each side of the equation, we obtain the following system: 

(3.62) 0=Q(F ,F ), 

(3.63) 8 t F + P ■ d s F = Q(F , Fi) + Q(F U F ), 

d t F 1 + P ■ 8^ = Q(F , F 2 ) + Q(F 2 , F ) + Q(F 1 ,F- L ), 



d t F 5 + P ■ d m F 5 = Q(F , F e ) + Q(F 6 , F ) 



E 

i+j=6 
l<i, j<6 



Q(Fi,Fj 



while the remainder satisfies the equation 
(3.64) d t F R . e + P ■ d s F R , £ - ~{Q(F ,F R . e 



6 

E 

i=l 



Q(FB,,e,F )}=e 2 Q(F R , e ,F R , e ) 
e 1 - 1 {Q{F il F R . e ) + Q(F R . e , Fi)} + e 2 A, 



with 



A d = l -e{d t F e + P-d s F 6 }+ £i+3 ~ 6 Q(F l ,Fj). 

i+j>6 
l<i,j<6 

By (1.19), equation ( |3.62[ ) implies that Fq must be a relativistic local Maxwellian 
Fo(t,x,P) = M = M( n(t, x ), 0(t, x),u(t, x);P), as in (1.17). Consequently, the 
remaining equations in (3.63) and below involve the linear operator: 



L(h)^-{Q(h,M) + Q(M,h)}. 

We remark that L is an integral operator involving only the momentum space 
variables. Furthermore, L is a linear Fredholm operator that can be inverted as 



long as the inhomogeneity (i.e., the terms in (3.63), and the equations below it, 



which are not of the form L(Fi)) is perpendicular to the five dimensional null 
space of the adjoint operator : Null(L^) = span{0i, . . . , ffc} = span{l, P, P }. 
This null space can be seen easily from the standard pre-post change of variables. 
In the preceding discussion, the notion of perpendicular and adjoint is the one 



corresponding to the usual L 2 momentum space inner product defined in (1.5). 
The operator L has the null space Null(L) = span{.A/f, MP, MP }. 

The aforementioned perpendicularity conditions can be checked by direct calcu- 
ion, 

(3.65) 



lation, so that equation (3.63) and the one below it imply that 
^i,d t F + P ■ d s Fo)_ = 0, (i = 1, • • • , 5), 

F x = -L- 1 (d t F a + P ■ 8 S F ) + $1, 



(3.66) 



1, 



^ i ,d t F 1 +P-d s F 1 -Q(F 1 ,F 1 ))_=0, (i 

F 2 = -L- 1 (d t F x + P • M - Q{F U F{)} + $ 2 . 
Above $1 and $2 are elements of the null space of L (i.e., L($i) 



,5), 



£(* 2 ) - 0). 

Applying the operator <9„ to each side of (jl.20[) , differentiating under the integral, 



and using (3.10d), it follows that equation (3.65) implies that the relativistic Eulcr 
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equations (1.22) are verified by the macroscopic quantities nLVf], 0[yW], and u[M] 
corresponding to the Maxwellian M = Fq. This explains the fact that in order to 
initiate the Hilbert expansion, we solve (with the help of Theorem^ for a smooth 
solution to the rE system using the variables (n,9,u). 

As discussed in Cercignani-Kremer [9j Section 5.5], the parameters in the expan- 
sion of $i in terms of the basis {M, MP , MP 1 , MP 2 , MP 3 } satisfy a linearized 
inhomogeneous version of the relativistic Euler equations; with the help of the ex- 



pression for Fi in (3.65), this enables us to find a solution F\ to equation (3.63). 
Furthermore, the higher order correction terms F%, F3, . . . , Fq can be solved for in 
the same way, where the corresponding inhomogeneous terms in the linearized rela- 
tivistic Euler equations depend upon the previous terms in the expansion. We refer 
to Cercignani-Kremer |9| Section 5.5] for more details on these terms in the expan- 
sion. We remark that a careful treatment of the non-relativistic Hilbert expansion 
is found in [8 24 25j. In particular, we are using the argument from |8j. These 
arguments carry over directly (once one identifies the null spaces in the relativistic 
case, as we have done above). 

We will use the following results, which are not studied in detail here. First, 
the terms F\, . . . , F§ are smooth in (t,x), and they also have decay in the mo- 
mentum variables, P. Consider, for example, F\. For the Newtonian version of 
and Caflisch [8] argue that L~ 



Grad 



25 



1 preserves decay in the momentum 
Their argument carries over directly to our case of the relativistic Boltzmann 
equation as follows. We can combine the argument in [8 25 with the relativis- 



tic estimates \ 1.15 1 

T 



Lemma 3.13 Lemma 3.14 and arguments as in 51 



to see 



that indeed L~ preserves momentum decay. This leads to the conclusion that 
— (jJtFo + P ■ dxFo^j decays at infinity as fast as M q for any < q < 1. Thus, 

Fx will decay as fast a,s M q , and F\ is smooth in (t, x) since the parameters in the 
expansion of $1 in terms of M, MP , MP 1 , MP 2 , MP 3 solve linear equations 
with forcing terms coming from the smooth functions n, 9, u. This argument is 
similar for the higher order terms in the expansion. 

The next proposition summarizes the estimates that we use in the next section. 

Proposition 3.9. Let (n(t,x),9(t,x),u(t,x)) be a smooth solution (see Remark 
1.8) of the rE equations (1.22) on a time interval [0, T] x M. 3 . Form the relativistic 
Maxwellian Fq — M(n,9,u; P) as in (1.17). Then the terms Fi,...,Fq of the 
Hilbert expansion are smooth in (t,x) € [0,T] x R| and for any < q < 1, they 
have momentum decay given by 

\F J (t,x,P)\<C(q)M q (n(t,x),9(t,x),u(t,x);P), (j = l,2,--- ,6). 

The constants in this bound are independent of (t,x, P). 

3.8. Relativistic Boltzmann estimates. In this section we prove our main re- 
sult, Theorem [2] Using Theorem [TJ we may assume that there is a sufficiently 
smooth solution (n, 9, u) to the relativistic Euler equations satisfying all of the 
desired properties in Theorem [2j We can then construct the local relativistic 
Maxwellian M(n(t,x),9(t,x),u(t,x); P) as in (1.17). After the analysis of Section 



3.7[ t he main point left is to estimate solutions to the equation for the remainder 
(3.64). We will outline the main strategy for these estimates after the statements 



of Lemma 13.101 and Lemma 13.111 below. 
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It will be useful to express the remainder as 

(3.67) = F R;e /VM . 

We use the notation ||/||2 = ||/||z / 2 (r3 xm 3 -) throughout this section. We define the 
linearized relativistic Boltzmann collision operator around M by 

C{h) = -M~ 1/2 {Q (VMh, M\+Q,(m, vCw/i) } . 

We also define a nonlinear operator by 

(3.68) r(/i,/) = M~ l ' 2 Q(yMh,VMf) . 

We recall the notation from Section ITTTl 

We further define the weighed L 2 (M. 3 i x M^) "dissipation" norm by 



\h\\l ^ 



dx / dP v(P) \h(x,P)f 



Above the the "collision frequency" , v(P) = v{J){P), is given by ( 1.15 ). Recall the 
weight function (1.6). We will sometimes write w = w(P) == wi(P). Furthermore, 

(3.69) h F ''" 



F R , s {t,x,P)/Jj(P) 



It will then be sufficient to estimate H/^^W and ||/i e ||oo to conclude Theorem 
[2] We prove the needed estimates in Lemma [3. 10| and Lemma [3. 11| just below. 

Let P denote the orthogonal L 2 (M. 3 p) projection with respect to the null space 
of the linear operator C, which is 

[VM, P 1 \/A4, P 2 Vm, P 3 Vm, p°Vm} . 

We know from e.g. [18 20 45] that there exists a number So > such that 
(3.70) (£h,h)p>S \\{I-P}h\\l. 
We will furthermore use the following L 2 - L°° estimates. 



Lemma 3.10. (L Estimate): We consider a smooth solution (see Remark 1.8) 
(n(t,x),6(t,x),u(t,x)) to the relativistic Euler equations (1.22) generated by Th 



eo- 



i[7} LetM(n,9,u;P), f e , h s be defined in 



and let 8q > be as in the coercivity estimate (3.70 
e > and C = C{M,F Q ,F X , 



3.67), and (3.69) respectively, 
). Then there exist constants 



, F 6 ) > 0, such that for all e € (0, Eq) we have 

|ll/ £ |ll(*) + |||{I - P}f\\l(t) < C{V~e\\e 3/2 h F \U e (t) + 1} {11/111 + ll/ E || 2 } ■ 

Above and below the constant C(A4, Fq, Fx, . . . , Fg) depends upon the L 2 norms 
and the L°° norms of the terms M., F , Fx, . ■ . , F 6 as well as their first derivatives. 



Lemma 3.11. (L°° Estimate): Under the assumptions of Lemma 3.10 there exists 
£0 > and a positive constant C = C(M, Fo, Fx, ■ ■ . , Fq) > 0, such that for all 
£ £ (0, £q) and for any £ > 9 we have 



sup ||e 

0<s<T 



3/2 ^||oo,£( S ) <c(||£ 3 / 2 Moo/+ SUp ||.f || 2 ( S )+£ 7 /2 
L 0<s<T 
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As we will soon explain, these two lemmas together imply Theorem [5] These 
estimates are motivated by the L 2 — L°° framework from |30| . Similar lemmas have 



also been used to study the non-relativistic Hilbert expansion in 31 . Indeed, the 
short proof of Theorem [2] is essentially extracted from [31] , modulo these lemmas. 

The main strategy is as follows. We first control the remainder equation (3.64) 
as in Lemma 3.10 using the L 2 energy estimates from 



45 



such as Lemma 3.12 



below. To finish the proof of Theorem [2] we also need L°° estimates such as those 
in Lemma 3.11 These can be proven using Duhamel's principle (3.74 1 and further 
estimates from [45] such as Lemmas |3.13| [3~14 and 3.15 below. In this framework, 
the key idea is to control the solution in L°° by the L 2 norms of the solution and 
the L°° norm of the initial data. 



Proof of Theorem [J| We will use the main estimates in Lemma |3.10| and Lemma 
|3.11| After applying the standard Gronwall inequality to the differential inequality 
in Lemma 3.10 we obtain 



where 



iirii^) + i<c(nri|2(o) + i) e cta ^\ 



a(T,e) = ^e sup He 3 ^! «,,«(*) + 1- 

0<t<T 



By Lemma [3TT] a(T, e) < b(T, e) with 

b(T,e) = ^(\\e 3/2 h Q \U e + sup ||/ s (s)|| 2 + e 1 ' 2 

\ 0<s<T 

We have shown ||/ e || 2 (t) < C (||/ || 2 + 1) e ctb(T ' e) . Notice that 

e x <C(l + x), if 0<a;<l. 
Therefore, for e sufficiently small, on some short time interval, we have 

||/ e || 2 (<)<C(||/o||2 + l)(l + v / ifl|£ 3/2 /io||oo/+ sup ||/ £ ( S )|| 2 

I \ 0<s<T 

Hence, there exists £o > such that for < e < Eq we may conclude 
sup |!/ £ (s)|| 2 < C T {1 + H/0II2 + ||£ 3/2 Moo,a. 

0<s<T 

This procedure works for some short time interval, and then the inequality above 
follows in general by a continuity argument. This last estimate and Lemma |3.11| 
together imply the main estimate in Theorem [2j □ 



For the remainder of this paper, we will discuss the proofs of the Lemmas 3.10 
and |3~TT) To this end, we will use the following nonlinear estimate. 

Lemma 3.12. For any £ > 9, we have the following estimate for the collision 



operator (3.68): 

(3.71) \(T(h 1 ,h 2 ),h 3 ) P \ < CH^Hoo/IIMHMb- 
Furthermore, if \ is an V rapidly decaying function, then we have 

(3.72) \(r(h 1 ,x),h 3 )p\ + \(r(x,h 1 ),h 3 ) P \<C\\h 3 \\4h 1 \\ v . 
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In the proof of this lemma below, we only require that \ satisfies the rapid decay 
condition |x(P)| < C(P°)~ m , with m > 3. On the other hand, in our applications 
below, we will consider smooth functions x with exponential decay, which is a prop- 
erty possessed by the relativistic Maxwellians defined in (1.17). 



Proof of Lemma 3.12 We notice from (|3.68|) and (|1.39|) that 

|r(/n,MI < I dQ 



doj v a{ Ql d) e' aQ ° \hi(P')h 2 (Q')\ 
+ [ dQ [ du v o-(q,#) e~ aQ ° \h x {P)h 2 {Q)\ =1 + 11. 

JR3 JS 2 

We estimate first the piece without post-collisional velocities, denoted 77 above: 



\(H,h 3 )p\< 



dP 



dQ 



< \\h 



3 Woo.e 



dP 



dQ 



dcu v o-(q,#) e~ aQ ° |/ii(P)/i 2 (Q)/i3(P)| 
duj v a(g, $) 



S 2 



|/ti(P)/t 2 (Q)| 

(Q a y(p°y ' 

Above, we made use of the trivial estimate e~ a Q° < (Q°)~ i . Applying the Cauchy- 
Schwarz inequality, we conclude that 

1/2 



\(n,h 3 )p\<\\h 3 \\ OOJ n 

i=l,2 



dP 



dQ 



/ \ \hi(P)\ 

°(e>*) iQ oy { poy 



We use the (P, Q) symmetry to interchange the values of P and Q in the integral 
involving h 2 above. Since t > 9, we observe from our hypothesis above (1.15) that 

v a{g,fl) 



(P U ) 



o\-e 



dQ 



duj 



(Q°Y 



< 1. 



Estimates of this type are proven for instance in |45[ Lemma 3.1]. ^From here, the 
first estimate in Lemma |3.12| follows for term 17. Similarly, for term / we have 



3/P 



< 



dP 



dQ / duo v a(g, d) 



s 2 



< H^siu,! n 

i=l,2 



dP 



dQ I du> v a(g, $) 



\hi{P')h 2 (Q')h 3 {P)\ 

- ,2 \ 1/2 



(Q°Y(p°Y 



Above, we used the (P',Q') symmetry to interchange the values of P' and Q' in 
the integral involving h 2 (P') above. By the pre-post collisional change of variables, 
which is dPdQ = p0 qo dP dQ , we have 



dP 



dQ / duj v a(g, $) 



MP') 



s 2 



(Q°Y(p°Y 



= / dp 



dQ 



, f a, H P )\ 



Above, we used the fact that that the kernel of the integral is invariant with respect 



to the relativistic pre-post collisional change of variables from 19 
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Now it can be seen that P° < Q' Q P'° and Q° < Q'°P'°. This is the content of 



20 Lemma 2.2]. From here, since I > 9 we observe that there exists a small 5 > 



such that 



duJ v t a ±i*i < i. 



(Q 



o\e/2+s 



This establis hes (|3.71 |. 

To prove (3.72), we choose a cut-off function \ satisfying x(P) < (p(^-m f or 
any m > 3. We will estimate the 7 and 77 terms from the top of this proof with 
hi = x m th" 3 nrs t case. In particular, as before, we have that 

\(n,h 3 )p\<C [ dP [ dQ [ duv t a{ e ,<f) e- aQ ° Ih^hiP)] 

JR3 JR3 Js 2 



We also have 



\{I,h 3 ) P \<C [ dP f dQ [ duv a(Q,§)e-^°{Q'°)- m \h 1 (P')h 3 (P)\ 

JR3 JR3 JS 2 

<c n (7 di 

,•£13 



dQ / dw V0 <7( e ,tf) (g°)- m (Q'°)- m ^(P') 



where we have used similar reasoning as we did for the previous terms. After the 
pre-post change of variables, as before, we conclude that h 3 )p\ < C||/i3||„||/ii||„. 
as desired. We have thus shown that |(r(/ii, x)> h 3 )p\ < ||/i3||i/||/ii||i>- The last esti- 
mate involving the term T(x, hi) follows in exactly the same way. □ 

With this lemma, we are ready for the 



Proof of Lemma 3.10 Since Fp. e satisfies (3.641, the function f £ from (3.671 satisfies 
the equation 



£ \ VM 



^ e*- 1 {r (Fi/VM, f) + r (> , f/Vm) } + e 2 A 



i=i 



where 



A = -eM- 1/2 {d t + P ■ d s }F 6 + ^ £ i+j ~ 6 T (fJVM, Fj/VM 



i+j>e 

l<i,i<6 



We now take the L 2 inner product (in Kj x Rp) of this equation with f e . It follows 
from (I3.70D that 



'hi- + r ■<>,/- + - £ c(n,r) _ > I|||/-||2 + ||| { i_p} /£ 
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In the remainder of the proof we will give upper bounds for the other terms. We 
point out that M~ 1 / 2 {dt+P-d x ^^fM is a first order polynomial in p. Furthermore, 



{d t + p ■ d s }VM 
Vm 



< Wl (P) C(n,u,0). 



Here C — C(n, u, 9) depends upon the the L°° norm of the first derivatives of 
the fluid variables. Choose a = f 2/(4 — b), where b is any number subject to the 



constraints described just above (1.15). With that, for any k > 0, we obtain 



{dt + P- d s }VM 
\/M 



r,r ) = I dxdp+ 

I x-P Jw ^ p )>7S J Wl (P)<- 

fe||2 



dx dP 



< C 2 (n,u,6)\\w 1 f e \\ L < x , iw{ p ) >j k) \\f e \\ L 2 + C 00 (n,u,8)\\ y /mf e \\ La(w(P) <K [ y 



Above, C2(n, u, 6) denotes a constant which depends upon L 2 norms of the first 
derivatives of the fluid variables and the L°° norms of the fluid variables alone. This 
estimate follows from the Cauchy-Schwarz inequality. The constant C^in, u, 9) also 
depends only on the L°° norms of the fluid variables and their first-order derivatives. 



Furthermore, from (1.39), (3.69) and (3.67), we observe that 



\Wx{P)\\f(P)\ <Wi- t (P)\ We (P)h £ \ < -^\ Wi (P)h S \, if Wl (P)>-r. 



The above estimate utilizes ^-l>8>(4-6) = 2/a and J < CM in (L39j) 
Additionally, since A = 2 — | , we have the following weight estimates: 



,l/5 ffl V2-l<^!^/2-l<^(pO)V2 if w (p)<JL. 

~ e e v ; v ' ~ e a 



With these two computations, we estimate 



{dt + P- d x }VM 



< 



-^jl II h e \\ oo.t || f || 2 + Coo (n, u, 9) || vW e II mw(P)< % ) 



< C K s 2 \\h e \ 



■ c\\v^pn\i Hw{P) ^ } 



<C K£ 2 ||/ l £ || 00 ,,||/ £ || 2 + C'||/ £ || + c'^||{I-P}/ e ||2. 



Choosing n to be sufficiently small, these are the estimates that we will use for the 
term involving derivatives of the local Maxwellian M. 



We use Lemma 3.12 and (3.69) to obtain for I > 9 that 



Hnrj £ )j e ) S ;p\ <y/i\\e 3/2 h*\u t \\n 
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Above we have also used (1.39) with J < CM. We utilize the second estimate in 
Lemma [3 . 1 2 1 and Proposition |3.9| to achieve 

t ^ { (r {Fi/VM, f) , f ) s: _ + (r (f, Fi /VM) , r)\ 



< IIP/ 6 



II{i-p}/T<II/1| 2 + -II{i-p}/X- 



We have used Proposition 3.9 (1.16), and (1.39) to conclude that Fi/vM. is a 
rapidly decaying function as in the statement of Lemma |3.12| Similar to the above 
estimates, we have (e 2 A,f £ )p < £ 2 ||/ £ ||l 2 ^ II/ £ ||l 2 - Note that | = 1 was added 
above so that we can obtain the estimate in Lemma |3.10| by absorbing this term into 
the l.h.s. of the inequality for e > sufficiently small. In particular, we conclude 
our estimate by choosing k small and then supposing that e < J^. □ 

We are now ready to consider the L°° estimate for h £ in Lemma [3.11| We first 
expand 

- J- 1/2 {Q(M, VJh) + Q(VJh, M)} = v{M)h - K(h), 

where K = K 2 — K\. This will be an important term in the equation for the 
remainder (3.64) once we plug in the ansatz (3.69). This is observed in the proof 
of Lemma 3.11 The operators Ki(h) and K 2 (h) are defined as 

Kx{h) = J- 1/2 Q-{M,VJh), 

K 2 (h) = J- 1/2 {Q + {M,VJh) + Q+(VJh,M)}, 

while 

u(M) = J~ 1/2 Q~(VJ,M) = Q~(1,M). 

Above, the operators are the usual gain and loss parts of Q from (3.4). More 
specifically, the operators Ki can be expressed as 



K^h) 



K 2 {h) 



R 3 xS 2 



duidQ v a(g, ■&) 



dudQ v a{ Ql d) <M(Q') 



dudQ v a(g,d) \M{P') 



sfW) 



MP') I 
h(Q') 



Given any small number 77 > 0, we choose a smooth cut-off function x — x(q) 
satisfying 

1 if q > 277, 
if q < r/. 



xio) 
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This Lorentz invariant cut-off function was previously used in |45| . We now define 



(3.73) 



(h) d - f / dwdQ (l-x(e)) v a{Q,$) M(Q')^3^h{P>) 
dujdQ (l-x(e)) v, v(q,0) M(P')^Qh(Q>). 



We also define K\ x {h) in the same way. We will use the splitting K = K 1 X +K x . 
We will use the Hilbert-Schmidt form for K x , it is given by 



K x (h) 



dQ k x {P,Q) h(Q), (i = 1,2). 



This form is computed explicitly in [45j Appendix]; see also [12]. Since the closed 
form expression for the kernel k x (P, Q) can be quite complicated, we simply state 
the following useful estimate: 

Lemma 3.13. 45] Lemma 3.2]. There exists a constant Q > such that the kernel 
enjoys the estimate 



k x (P,Q)\ < (P°Q°y C (d + (P° + Q Q )- (P) - /2 )e- cl 



■\P-Q\ 



where (/?)- = max{— /3,0}, (3 is the parameter, and C\ is the constant from above 

Note that in the notation of this paper we use P° to denote the quantity which is 
called po in [45J . The estimate above is proved in [45[ Lemma 3.4] for the case of soft 
potentials, i.e. (/?)_ = b, a = 0, and C\ = in ( 1.15 ). However the generalization to 
the case above is immediate and follows directly from the proof. In the case above, 
we have the exact estimate ( = min{l — (a + (j)-)/2, min {2 — |7|,4 — 6, 2} /4}. 
We will now quote the estimate of the operator K l ~ x as follows: 

Lemma 3.14. [45[ Lemma 4.6]. Fix £ > 0. Then given any small rj > 0, we have 



\w t {P)K^{h{P))\<y 1 e-^ p °\\h\\ 00 . 

Above the constant 7 > is independent of rj. 

We will also use the following nonlinear estimate: 

Lemma 3.15. [45[ Lemma 5.2]. For any I > 0, we have the following L°° estimate 
for the nonlinear Boltzmann collision operator: 

w e (P)J- 1/2 (P)Q (hxVJ, h 2 VJ^) (P) 

We are now ready to prove Lemma |3.11| 



< 



v{P)\\h x 



l 2\\oo,l 



Proof of Lemma 3.11\ From (|3.64|) and (|3.69|), we obtain 
d t h e + P-d s h 



v(M) 



'J 



6 

i=l v J 
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with 



A d A c -Z{d t + p. ds }F 6 
v J 



E 



.j+i-6 



i+j>6,i<6,j<6 



<J 



We will prove Lemma 3.11 by iterating twice the representation of solutions to this 
equation in terms of the Duhamel formula. 

For purposes of the proof, we define yi == X — P(t — s). Let K = K x ~ x + K x be 
the splitting from (3.731 and define 



v(t,s) 



dr v{M){r,x) 



Furthermore, using (1.39) we have v(J)(P) ps v(M){P). ,/From this and (1.151, we 
have that 

u(t,s)^y(J)(P)(t-s)^(P°f/ 2 (t-s). 
Here and below A w B indicates that 3 C > 1 such that < B < CA. 



Now by the Duhamel formula we have 

j7(t,0) 



(3.74) h e (t,x,P) = cxp - 



ds exp 



h £ (x-Pt,P) 
v(t, s 



(-T 2 ) 



ds exp 



^ 1 -^(/i £ )( Sl2/1 ,p) 
^(^)( s , yi ,p) 



,/. s ( -^4=C(^VJ,/iVJ) (s,yi,P) 



as exp ^ 



ds exp(-^M)e 2 i( s , 2/1 , J P). 



We will now simply write s) as v{t — s). We will use the following basic estimate 
several times below 



ds exp I 



< 



We will now estimate each of the terms in (3.74 1. 

Given r\ > 0, we r ecall the splitting K = K r ^ L + K x as in p.73|. ^From Lemma 



3.14 



for any r\ > 0, the K 1 x term in (3.74) multiplied by wg is bounded as 

' A I W \ 

ds exp | 

y(t-s)' 



< 



ds <^(- l ^^')K 1 - x (h e )(8,y 1 ,P) 
- sup ||ft. £ (t)||oo.<! / ds exp ( 

^ 0<*<T Jo V 



-)v <<ne sup 11^(4)11 

0<t<T 



Next, we use Lemma [3.15| to conclude that 



we 



7i e VJ, /iVJ 



< 



"(-P) ll^ll^ 
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Then the fourth line in (3.741 is bounded by 

(3.75) Ce 2 f ds expf- \\h e (s)^ e < e 3 sup \\h E (s 
Jo \ e J o<s<t 

For the following terms in (3.74), from Lemma 3.15 again, we have 

6 



-i we 



- {Q (F u h £ Vj) + Q (h e y/J, Fi) } 
<cf^ v{P) \\h 



Fi/VJ < Cv{P) \\h £ 



\oo,e- 



We have used Proposition 3.9 and the upper bound in (|1.39|) to conclude that 



oo.i 



< C. Thus, the fifth line in (3.74) is bounded by 



(3.76) 



ds exp 



v{t - s) 

e 



0<t<T 



As above, the last line in (3.74) is clearly bounded by Ce . 

We have now estimated all the terms in (3.74) save one, which we denote by J 7 , 
and which is defined by 

ft 



T dcf 1 

7 ~ e 



Collecting the above estimates, we have established 



sup ||e 3/2 /i e W||oo,£ < Ce sup We^h^t^^+Ce 3 sup ||e 3/2 /i e e +Ce 3 

0<t<T 0<t<T 0<t<T 



+ C\\e 3 / 2 h \\ oo , e + Cw e (P)e 3 / 2 I 7 



To bound 7 7 , we will input h E in the form of (3.74) back into J 7 just below. 

We recall that the kernel of K x (h) is k x (P,Q). With this notation, it follows 
that 

L t = I ds expf- -i-i-U f dQ k x (P,Q) h e (s, yi ,Q). 
Jo v £ ' £ Jr 3 

We plug (3.74) for h e into the above to obtain 
rt 



(3.77) L ( = I ds exp 



o 



v(t, s)\l 



x exp 



e / e 



dQ k x (P,Q) 
hl(x - P(t - s) - Qs, Q) + 7 7)7 + H-y. 



We also introduce the notation y2 = x — P(t — s) — Q(s — s'). Furthermore, we 
define 



J 7 . 7 = f ds exp ( - KP)(M) ) - f dQ fc*(P, Q) 

Jo V £ ' £ Jm 3 



x / ds' exp ( — 



K* (h s )(s',y 2 ,Q). 
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Additionally, the term H 1 can be expressed as 



(3.78) ff 7 



J ds exp ^ 



ds exp I 
ds' exp 

i7(P)(M)\l 



i7(P)(i, a )\l 



c?s exp f 

6 



7J 



e / e 
/ k x (P,Q) f ds' exp( 

Q(V\/j, ft £ Vj) (s',2/ 2 ,Q) 



?'(g)( S;S ')' 



dQ k x (P,Q) j ds' exp( 

E £i_1 Jj {Q (Fi,^) + Q (h E VJ, Ft) } (s', y 2 , Q) 
V{P){t,s)\l 



v'{Q){s,s'Y 



i=l 



ds exp 



( 



£ ./»3 



dQ k x (P,Q) 



J ds' exp ^ 



?'(Q)(s,s') 



e 2 A(s',y 2 ,Q). 



Above we are using the necessary additional notation 



v'{Q){s,s') 



dr v(M){t,x). 



We have v(J)(Q) s» j/(yVl)(<2), which implies 

^(Q)( a y)«i/(j)(0)( a - a , ) 1 

which we will now simply write as v(s — s') as we did in the previous case of v. 

Using similar arguments to (3.75) and (3.76), we can control all the terms in 
(3.77) and (3.78) except the second term in (3.77) by the following upper bound: 



C{ WKW^t + e 6 sup Wh^sW^+e sup Wh^s)^ + Cs< 

0<s<T 0<s<T 



We now concentrate on the second term in (3.77), I 1 



7.7- 



We claim that for any small rj > 0, the following estimate holds: 

\w t I^\<T)' SUp H^OOHoo^ + Cy SUp ||/ £ (s)|| 2 . 
0<s<T 0<s<T 



This is proved in 45 Lemma 4.4]. We note that |45l Lemma 4.4] is explained in 
detail for the soft potentials and for general momentum weights parametrized by 
k > 0. Here we only use the k = case from 45 Lemma 4.4]. Furthermore, this 
estimate can be easily extended to the full range of hard and soft-potentials. The 
proof for the hard potentials follows in exactly the same way as the estimate for 
the soft potentials, but several technical simplifications of the proof are possible in 
the hard potential case. 
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With this last estimate we will now finish the proof. We first collect all of our 
estimates as follows: 

sup \\e 3 / 2 h%s)\\ 0O . i <C{r 1 + r 1 '} sup \\eV 2 tf{s)\\ oo4 + e 7 ' 2 C 

0<s<T 0<s<T 

+ CJe 3 / 2 h \U t + s^ 2 C sup \\e 3 / 2 h% s )\\l e + C n , sup ||/ £ (s)|| 2 . 

0<s<T 0<s<T 

We now choose rj and rj' small in such a way that C{r]+rj'} < \ . Then for sufficiently 
small e > 0, we obtain 

sup \\e i ' 2 h s { s )\\ oo4 <c{\\e^ 2 h \\ ool + sup ||/%s)|| 2 + e 7 ' 2 \ , 

0<s<T I 0<s<T J 

and we conclude our proof. □ 
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